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Abstract 

It is by now well-known that one can recover a potential in the wave equation from the knowl- 
edge of the initial waves, the boundary data and the flux on a part of the boundary satisfying 
the Gamma-conditions of J.-L. Lions. We are interested in proving that trying to fit the discrete 
fluxes, given by discrete approximations of the wave equation, with the continuous one, one re- 
covers, at the limit, the potential of the continuous model. In order to do that, we shall develop a 
Lax- type argument, usually used for convergence results of numerical schemes, which states that 
consistency and uniform stability imply convergence. In our case, the most difficult part of the 
analysis is the one corresponding to the uniform stability, that we shall prove using new uniform 
discrete Carleman estimates, where uniform means with respect to the discretization parameter. 
We shall then deduce a convergence result for the discrete inverse problems. Our analysis will be 
restricted to the I-d case for space semi-discrete wave equations discretized on a uniform mesh 
using a finite differences approach. 

Key words: Inverse problem, Discrete wave equation, Discrete Carleman estimate, Stability, 
Convergence. 

AMS subject classifications: 35R30, 35L05, 65M32, 65M06 

1 Introduction 



In this article, our goal is to study the convergence of an inverse problem for the f-d wave equation. 
Before introducing that problem, we shall present which inverse problem we are dealing with in 
the continuous setting. 
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The continuous inverse problem. For T > 0, we consider the following continuous wave 
equation: 

Buy - d xx y + qy = g, (t, x) e (0, T) x (0, 1), 

y(t,0)=g°(t), y(t,l)=g 1 (t), te(0,T), (1.1) 
2/(0,0=2/°, 9*2/(0, •)=W 1 . 

Here, y = y{t,x) is the amplitude of the waves, (y ,?/ 1 ) is the initial datum, g = q(x) is a 
potential function, g is a distributed source term and (g°,g 1 ) are boundary source terms. 

Of course, this problem is well-posed in some functional spaces, for instance: If (y ^ 1 ) G 
ff^Ojl) x L 2 (0, 1), 3 G L'fO.TjL^O, 1)), g i G fl^O.T) for i = 1,2, with the compatibility 
conditions y°(0) = ff°(0) and y°(l) = g\Q) and q G L°°(0, 1), the solution y of (fl~Tj) belongs to 
C([0, T]; -ff 1 (0, l))nC 1 ([0, T], I/ 2 (0, 1)). Such result is well-known except perhaps for the condition 
on the boundary data, which is a consequence of a hidden regularity result and a duality argument, 
giving a solution of (|l.ip in the sense of transposition - see detailed for instance in [27] . Under 
this class of regularity, using again a hidden regularity result in [2S], we can prove that d x y(t, 1) 
belongs to L 2 (0,T). 

We can therefore ask if, given (y°, y 1 ), g, (g°, g 1 ), the knowledge of the additional information 
d x y(t, 1) for a certain amount of time allows to characterize the potential q. We emphasize here 
that the data (y ,?/ 1 ), g, (y ,'? 1 ) are supposed to be known a priori. 

It has been proved in pQ that this question has a positive answer provided that T is large 
enough (T > 1 here) and y G -ff 1 (0, T; L°°(0, 1)). Of course, to guarantee this regularity without 
any knowledge on q, we may impose some stronger conditions on the data (y°, y 1 ), g, (g°, g 1 ), see 
e.g. in Remark 11.21 below. 

Let us precisely recall the results in pQ. For m > 0, we introduce the set 
Lf m (0, l) = {qe L°°(0, 1), s.t. \\q\\ L oo (0A) < m}. 

It will also be convenient to denote by y[q] the solution y of (|1.1[) with potential q. Assuming that 
p G L^? m (0, 1) is a given potential, we are concerned with the stability of the map q i-> d x y[q](-, 1) 
around p. Then we have the following local Lipschitz stability result: 

Theorem 1.1 ( F). Let m > 0, K > 0, r > and T > 1. 

Let p in L^ ) m (0, 1). Assume that the corresponding solution y[p] of equation (|1.1[) is suc/i t/iat 

l|y[P]llijl(0,T;£oo {0il)) < (1.2) 

Assume also that the initial datum y° satisfies 

inf{|»°(x)|,ar6(0,l)}>r. (1.3) 

Then for all q G L5f m (0, 1), 9t x y[p]{-, 1) — dt x y[q](-, 1) G L 2 (0, T) and there exists a constant C > 
that depends only on the parameters (T,m, K,r) such that for all q G L5f m (0, 1), 

\\dt x y\p]{-,l) - dt x y[q]{-,l)\\ L 2 {0 , T) < C \\p - q\\ L 2 (0>1) , (1.4) 
ll9-PL=(o,i) < C\\d tx y\p}{-,1) -d tx y[q](;i)\\ L 2 i0tT) ■ (1.5) 

Estimate (|1.5p gives the Lispchitz stability of the inverse problem and (|1.4[1 states the con- 
tinuous dependance of the derivative of the flux of the solution with respect to the potential. 
Together, these two estimates indicate that the above result is sharp. Note however that estimate 
(|1.4p is, by far, the easiest one to obtain. 

Remark 1.2. The condition (|1.2[) can be guaranteed uniformly for p G L< m (0,l) with more 
constraints on the data (jAy 1 ), g, (y°,y 1 ) m (|1.1|1 . for instance: 

(V°,» 1 )6ff 3 (0,l)xff 1 (0,l) 1 

g G W /1,1 (0, T; L 2 (0, 1)), (ff ,;? 1 ) G (iJ 2 (0,T)) 2 , 
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under the compatibility conditions 

g°(0) = y°(0), g\0) = y°(l), d t g°(0) = ^(0) and d t g\0) = y\l). 

Indeed, under these assumptions, dty[p\ belongs to the space C°([0, T]; i? 1 (0, l))nC 1 ([0, T]; L 2 (0, 1)) 
(see JF7|9, with estimates depending only on m and the norms of(y°,y 1 ), g, (g ^ 1 ) in the above 
spaces. Therefore, due to Sobolev's imbedding, y[p] satisfies (|1.2[) for some constant K > that 
can be chosen uniformly with respect to p G L5f m (0, 1). 

The method of proof of Theorem 1 1.1 1 is based on a global Carleman estimate and is very close 
to the approach of [22], that concerns the wave equation with Neumann boundary condition and 
Dirichlet observation for the inverse problem of retrieving a potential. Actually, it also closely 
follows the approach of 33 but the work [1] requires less regularity conditions on y. 

The use of Carleman estimates to prove uniqueness in inverse problems was introduced in [9] 
by A. L. Bukhgeim and M. V. Klibanov. Concerning inverse problems for hyperbolic equations 
with a single observation, we can refer to [30] . [31] or [34], where the method relies on uniqueness 
results obtained by local Carleman estimates (see e.g. [20] . [26] ) and compactness-uniqueness 
arguments based on observability inequalities (see also [35]). Related references [25], [5T] and [23] 
use global Carleman estimates, but rather consider the case of interior or Dirichlet boundary data 
observation. Let us also mention the work [5j for logarithmic stability results when no geometric 
condition is fulfilled. 

Discrete inverse problems. In this paper, we would like to address the question of the 
numerical computation of an approximation of the potential p G L°°(0, 1), on which we assume 
the additional knowledge that its L°°(0, l)-norm is bounded by some constant m > 0. 

A natural approach is to find ph G L< m (0, 1), or rather in a discrete version of it denoted by 
L^ <m (0, 1) that will be made precise later (see (|2.6[) ). such that 

dtd x y h [ Ph ](t,l)~d t d x y[p](t,l), t€(0,T), (1.6) 

where y^ is the solution of a corresponding discrete wave equation with potential p^ (here, h > 
refers to a discretization parameter) and the meaning of (|1.6p has to be clarified. The question 
is then the following: Does ()1.6[) imply pu — p ? Or, to be more precise, can we guarantee the 
convergence of the discrete potentials ph toward the continuous one p ? 

Our analysis will focus on this precise convergence issue. To sum up in a very informal way 
our results, we will show that the convergence indeed holds true (Theorem 14. ip . provided a Ty- 
chonoff regularization process is introduced, and the key estimate is a stability estimate for the 
discrete inverse problem (Theorem l3.1[) . given by appropriate global discrete Carleman estimates 
(Corollary 12.41 and Lemma I2TT31) . 

To be more precise, for N 6 N, set h = 1/(N + 1), and let us consider the following semi- 
discrete 1-d wave equation: 

dttVj,h - {^hVhij + qj,hVj,h = 9j,h, t G (0,T), j G [1,7V], 

yo,h(t) = 9°(t), VN+i,h{t) = gl{t), te[o,T], (1.7) 

Vi,k(0) = Vlh, ftVJ,fc(0) = vlh, 3 G [1, N], 
where 

(A h y h ) 3 = ta (Ki+i.h ~ 2j/j,fc + Vj-i.h) 

denotes the classical finite-difference discretization of the Laplace operator and where {y^ hi v),h) 
are the initial sampled data at Xj = jh, g\ G L 2 (0,T), i = 0, 1 and G L 1 (0, T; L\{Q, 1)) are 
the boundary and source sampled data. Here and in the sequel, L\{Q, 1) denotes a discrete space 
endowed with a suitable discrete version of the L 2 (0, 1) -norm that will be introduced later (see 
(|2.4[) '). In the following, it will be important to sometimes underline the dependence of yn in (|1.7[) 
with respect to the potential qn- This will be done using the notation yh[qh]- 
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Note that, using these notations, the discrete normal derivatives of yh[qh] at the point x = 1 
is naturally approximated by 

yN+i,h[lh] - VN,h[lh] _ gl(t) yN,h[lh] 



h 



h 



h 



The fact that g\ is known will allow us to simplify the difference of the discrete normal derivatives 
of yh[qh] and yh[ph] simply as 

VN.h[<lh] _ VN,h[Ph] 

h h ' 

In order to prove the convergence of the inverse problem, we shall develop a Lax-type argument 
for the convergence of the numerical schemes that relies on: 

• Consistency: lip £ L< m (0, 1), there exists a sequence of discrete potentials ph £ L^ <m (0, 1) 
such that 



Ph — > P in L (0, 1), 

h— yO 

d t yN+i,h[ph\ — d t yN,h[Ph] 
h 



ftflUlM-,1) inL 2 (0,T). 



(1.8) 
(1.9) 



Uniform stability: There exists a constant C independent of h > such that for all 
(q h ,p h )&L%< m (0,l) 2 , 



\\q h -Ph\\ L 2 (01) 



< c 



dtyN,h[lh] d t yN,h[ph] 



li 



h 



(1.10) 



L 2 (0,T) 



Of course, the consistency is the easiest part of the argument and will be detailed in Section [4] 
The most difficult one comes from the stability estimate (|1.10[1 . 

Actually, as we shall explain below, we will not get (| 1 . 10|> . but we shall rather prove an 
estimate of the form 



hh -Ph\\ L 2 (0tl) < C 



dtVN,h[qk] d t y N ,h\ph] 



L 2 (0,T) 

+ Ch \\d£d t tyh[qh] - d^d tt yh[Ph] 



(1.11) 



where 



(d h Vh)j = 



Vj+i,h - y,,h 



for some C > independent of h > - see Theorem 13. II for precise statements. 

This is still compatible with the Lax argument: the added observation operator weakly con- 
verges to as h — > 0, since hd£ is of norm bounded by 2 on L\(0, 1), and obviously converge 
to zero for smooth data. Therefore, in the limit h — > 0, this term disappears and (|l.ll|l still 
yields ([T3j) . 

Of course, this should be taken into account into the consistency argument: Given p G 
L°°(0, 1), one should find a sequence ph such that (|1.8|l - (|1.9|) hold and 



hd^d tt yh[ph] 



in L 2 ((0,T) x (0, 1)). 



(1.12) 



We refer the reader to Theorem 14.21 for precise assumptions and statements concerning the con- 
sistency. 

The convergence result for the discrete inverse problems toward the continuous one is then 
given in Theorem 14.11 and takes into account the previous comments. Roughly speaking, we will 
prove that, given any p £ L°°(0, 1) and any sequence ph G L^< m (0, 1) such that the convergences 
(|1.9[) and (|1.12|) hold, the discrete potentials ph converge to p in L 2 (0, 1) as h — > 0. We refer to 
Section \4A\ for precise assumptions and statements. 
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The proof of the uniform stability estimate (|1 . 1 If) (see Theorem 13. 1|) is based on a discrete 
Carleman estimate for (II. 7[) . which should be proved uniformly with respect to h > (see Corol- 
lary [2]4]). This is the main difficulty in our work. 

First, a discrete version of the continuous Carleman estimate yielding the stability (|1.5f) cannot 
be true as it is. Indeed, that would contradict the results in 19, 25, 36, 16 that emphasize the 
lack of uniform observability of the discrete wave equations. This is due to the fact that the 
semidiscretization process that yields (|1.7[) creates spurious high-frequency solutions traveling at 
velocity of the order of h, see e.g. [321129] . Hence, they cannot be observed in finite time uniformly 
with respect to h > 0. 

We shall therefore develop a discrete Carleman estimate for the discrete wave equation (|1.7[) 
which holds uniformly with respect to the discretization parameter h > 0. We will use the same 
Carleman weights as in the continuous case. Though, the discrete integrations by parts will 
generate a term which cannot be handled directly. This will correspond to a term of the order 
of 1 at high-frequencies of the order of 1/h, whereas it is small for frequencies of order less than 
1/h, thus being completely compatible with the continuous Carleman estimates and the analysis 
of the observability properties of the discrete wave equation. One can see [361 116] for review 
articles concerning that fact. 

Uniform Carleman estimates for discrete equations have not been developed extensively so 
far. The only results we are aware of concern the elliptic case [6] [7j [8] for applications to the 
controllability of discrete parabolic equations, in particular in [8]. More recently in [14], discrete 
Carleman estimates have been derived for elliptic equations in order to prove uniform stability 
results for the discrete Calderon problems. 

Outline. The paper is organized as follows. Section [2] is devoted to the proof of discrete 
Carleman estimates for a 1-d semi-discrete wave operator. A uniform stability estimate for the 
related inverse problem is derived from it in Section [3] Convergence issues are finally detailed 
and proved in Section [4] and further comments are given in Section [5] 

2 Discrete Carleman estimates 

In this section, we establish uniform Carleman estimates for the semi-discrete wave operator. 

2.1 Continuous case 

We recall here the global Carleman estimates for the continuous wave operator. 

Let x° < 0, s > 0, A > and /3 G (0, 1). On [-T, T] x [0, 1], we define the weight functions 

ip = xp{t, x) and tp = ip(t, x) as 

^(t,x) = \x-x°\ 2 - fit 2 + C , ip{t,x) = e X ^ x) , (2.1) 

where Co > is such that tp > 1 on [-T,T] x [0, 1]. 

Let us begin by recalling the continuous Carleman estimate. This will make easier the com- 
parisons with the forthcoming discrete ones: 

Theorem 2.1 ( 1 ). Let Lw = d u w - d xx w, T > and f3 G (0, 1). 

There exist Ao > 0, so > and a constant M = M(so, Ao, T, /3, x°) > such that for all s > so, 
A > Ao and w satisfying 

Lw G L 2 ((-T,T) x (0,1)), 
w G L 2 (—T,T;Hq(0,1)), 
w(±T,-) = d t w(±T,-) = 0, 
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we have 



s\ / / ifie 2s<p (jdtw\ 2 + \d x w\ 2 ) dxdt + s 3 \ A I I ip 3 e 2s(p \w\ 2 dxdt 

< M f T f e 2sv \Lw\ 2 dxdt + MsX f <p{t, \) e 2av(tA) \d x w(t, 1)| 2 dt. (2.2) 

J-TJO J -T 

Carleman estimates for hyperbolic equations can be found in [20] and we refer the reader to 
the bibliography therein for extensive references. The Carleman estimate stated here can be seen 
as a more refined version of the one in 20, Theorem 1.2] but in the case of boundary observation 
and with the freedom on A. For the proof of this Carleman estimate, we therefore refer to pQ. 

Remark 2.2. Note that the above Carleman estimate holds without any condition on T. This 
might be surprising but this should not be since we assume that w(±T) = dtw(±T) — and 
therefore, the corresponding unique continuation result is: If w(±T) = dtw(±T) = 0, w 6 
L 2 (-T,T;H%(0, 1)), d u w - d xx w = and d x w(-, 1) = 0, then w = Q. 



2.2 Statement of the results 

In this section, we state uniform Carleman estimates for semi-discrete wave operators. 

Of course, since we work in a semi-discrete framework, the space variable x is now to be 
considered as taking only discrete values Xj = jh £ [0, 1] for j £ [0, . . . , N + 1] (recall that h = 
1/(N + 1)). Therefore, for continuous functions / (e.g. with ip, ip,...), we will write indifferently 
f(xj) or fj. We will also add the subscript h when we want to emphasize the dependence in the 
mesh size parameter h > 0, but we shall remove it as soon as the context clearly underlines that 
we are working for one particular h > 0. 

Moreover, by analogy with the continuous case, we will use the following notations: 

„ N , N - N+l 

/ fh = hY / hh, / fh = h^2f jlh , / fh = h^2f jth ,. (2.3) 

•'(0,1) j= i ■'[0,1) i=0 Ao,i] i= i 

Note that it also defines in a natural way a discrete version of the L p (0, l)-norms as follows: 
for p G [l,oo), we introduce L v h (0, 1) (respectively L?([0, 1))) the space of discrete functions fh 
defined for jh, j £ [1, NJ, (resp. j £ [0, iV]) endowed with the norms 

\\h\\l li0A) = [ \h\ P (^p. «, ([0il)) = / \h\ P ), (2-4) 
h -1(0,1) h J\o,i) 



and, for p = oo, 



IAIIl?°(0,1) = SU P \fi,h\ ( rcs P- II A II £~ ([0,1)) = SU P \fi,h\)- (2-5) 
h jeli.JV] je[o,iv] 



By analogy with L^ m (0, 1), we also define 

ih?<m(0, 1) = iqh = fe,fc)j6li,iV] G L^(0, 1), s.t. \\qh\\ L oo (0il) < mj . (2.6) 
We shall also use in the sequel the following notations: 

{m h v h )j = -i j i ; {mlv h ) 3 = {m h v h ) j+ i = J — ; 

la \ v j+l,h — Vj-lh /„+ \ /„_ > v j + l,h — v j,h 

(9hVh)j = 2h ; {dt v h )i = (d h v h ) J+ i = - J± - L j — — ! 

/ a „ s Vj+i,h - 2vj +Vj-i,h 

{/\ h vhh = j^p ■ 

One of the main results of this paper is the following discrete Carleman estimates: 
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Theorem 2.3. Let L h w h = d u w h - A h w h , T > and f) £ (0, 1). 

There exist so > 0, A > 0, s > 0, ho > and a constant M = M(so, A, T, e, /3) > independent 
of h > siic/i t/iai /or a/Z /i £ (0, /lo) and s £ (so,e/h), for all wh satisfying 

L h w h £ L 2 (~T,T;LU0,l)), 
w 0t h(t) = ui]v + i, h (t) = on (-T, T), 
w h (±T) = 9 tm (±T) = 0, 

uie /law 

s / / e 2sl/> |9 t TOft| 2 ctt + s / / e 2sv |9^"Wh| 2 + s 3 / / e 2sv \iVh\ 2 dt 

J —T J (0,1) i-T i[0,l) i-T i(0,l) 

<M [ T [ e 2 ^\L h w h \ 2 dt + Ms f e 28 ^ \(d^w h ) N+1 \ 2 dt (2.7) 

J -T J (0,1) J -T 

+Ms ( T f e 2sip \hd+d t w h \ 2 dt 

J -T J [0,1) 

for ip given by (|2.1|) . 

The proof of Theorem 12.31 will be given at the end of Section 12.61 
The following remarks are in order: 

• The weight function ip in the above discrete Carleman estimate is the same as for the 
continuous one. 

• In Theorem 12.31 the parameter A is fixed, whereas it is not in the continuous Carleman 
estimate of Theorem l2.ll Looking carefully at the proof of Theorem l2.3l one can prove that there 
exists A such that, for all A > Ao, there exist e(A) > and M = M(A) such that H2.7|) holds for 
all s > so (A) and sh < e(A). These dependences of e and M on A are very intricate and we did 
not manage to follow it precisely. 

• The fact that M is independent of h > is of major importance in the applications we have 
in mind. This is very similar to the observability properties of discrete wave equations for which 
one should prove observability results uniformly with respect to the discretization parameter(s), 
otherwise the discrete controls (obtained by duality from the discrete observability properties) 
may diverge, see e.g. [16] - 

• The range of s in Theorem 12.31 is limited to s < e/h. This is a technical assumption, that 
is not surprising when comparing it to [5] [7]. Indeed, for s of the order of 1/h, e S{p is a high- 
frequency function of frequency of the order of 1/h and therefore it does not reflect anymore the 
dynamics of the continuous wave operator. 

• A new term appears in the right hand side of H2.7|) . which cannot be absorbed by the terms 
of the left hand side. Though, this term is needed and cannot be removed. Otherwise, one could 
obtain a uniform observability result for the discrete wave equation, a fact which is well-known 
to be false according to [25]. Besides, this extra term is of the order of one for frequencies of the 
order of 1/h, whereas it can be absorbed by the left hand side for frequencies of smaller order. 
According to [16] . this indicates that the extra term in estimate (|2.7p has the right scale. 

Note that in the application we have in mind, we shall not use directly the Carleman estimate 
(|2.7p which involves the wave equation without a potential but rather one in which a L°° potential 
is allowed. Indeed, we have the following corollary: 

Corollary 2.4. Let T > and /3 £ (0,1). Let m > 0, q h £ L^< m (0, 1) and L h [q h ]w h = 
dttuth — A h Wh + qhWh- 

There exist so > 0, A > 0, s > 0, ho > and a constant M = M(so, A, T, m, e, j3) > such that 
for all h £ (0, ho) and for all s £ (so, e/h), for all Wh satisfying 

L h [q h ]w h £L 2 (-T,T;L 2 (0,1)), 
w>o,h(i) = w N+lyh (t) = on (—T,T), 
w h (±T) = 8 t w h (±T) = 0, 
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we have: 



if f e S{p \d t w h \ 2 dt + s f f e 2sv \d+w h \ 2 dt + s 3 f f e 2sip \w h \ 2 dt 

J — T J (0,1) J-TJ(0,1) J -T J (0,1) 

f e 2 ^\L h [q h ]w h \ 2 dt + Ms f T e 2av{t ^ \{dr,w h ) N+1 \ 2 dt (2.8) 

' J (0,1) J —T 



< M 



for ip given by (|2.ip . 



T J (0,1) 

+Ms f T f e 2sip \hd+d t w h \ 2 dt, 

J -T J [0,1) 



Proof. This is a simple consequence of Theorem 12.31 since Lhvjh = Lh[qh]wh — qn Wh, with qn 6 
L%< m (0, 1) leads to 



e 2s *\L h w h \ 2 dt<2 f f e 2s[p \L h [q h ]w h \ 2 dt + 2m 2 f f e 2s *\w h \ 2 dt. 

TJ(0,1) J -T J (0,1) J -T J (0,1) 



This last term can be absorb by the left hand side of (|2.7[) by choosing s large enough. This 
immediately yields (|2.8[) . □ 

Until the end of this section, we shall work for h > fix. We therefore omit the indexes h on 
the discrete functions to simplify notations. 

2.3 Basic discrete identities 

Below, we list several preliminary dentities that will be extensively used in the sequel. Let us 
begin with easy identities left to the reader: 

Lemma 2.5. The following identities hold: 

a 1 b 1 +a 2 b 2 ( eg + a 2 \ f &i + b 2 \ | h 2 / ai - a 2 \ f h - b 2 \ 



aibi — a 2 b 2 ( ai — a 2 \ ( b\ + b 2 \ ^ / ai + a 2 \ / &i - b 2 \ ^ ^ 



h \ h 

Using these identities, one can obtain the next lemma: 
Lemma 2.6. The following identities hold: 

h h 2 
ml=I + -dt\ m h = / + —A h = m^m h ; (2.11) 

dh = + dl ) = mtdl = d^ml = m^dl = d^ml ; (2.12) 

A, = 9+9" = d~dt ; (2.13) 

mt(uv) = (m+u)(m+v) + ^-(d+u)(d+v) ; (2.14) 

dfiuv) = (d£u)(mtv) + {m^u){dtv) ; (2.15) 
A h (pv) = (A hP ) (m h v) + 2(d h p) (d h v) + (m hP ) (A h v). (2.16) 

Proof. To begin with, one easily obtains (|2.1ip . since 



, + . Vj+i+Vj hv j+1 ~Vj h ,„ + , 

{mlv)j = J - = Vj + - -2— - — = V j + -(dtv)^ 



Vj+i + 2vj + Vj~i Vj+i — 2vj + Vj-i h 2 
(m h v) j = -2 ^ 3 — = Vj + -± -± 3 — =Vj + — (A h v) 

Similar computations left to the readers yield (12.121) and (|2. 13[) . 
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Identities (|2.14p - (|2.15|l are straightforward consequences of the formula of Lemma [2. 5 1 To get 
(|2.16|l . we do as follows: 

A h (pv) = 5^ (dt(pv)) 

= {(d+p)(m+v) + (m+p)(d+v)) 

= ( d h d h P)i m h m h v ) + ( m h d hP)( d h m h v ) + ( d h m hP)( m h d h v ) + ( m h m hP)( d h d h v ) 
= (A h p)(m h v) + 2(d h p)(d h v) + (m h p)(A h v). 

Note that this should of course be compared to the corresponding classical Leibniz formula 
A(pv) = vAp + 2Vp- Vt) + pAv. □ 

We now explain how discrete integrations by parts work: 

Lemma 2.7 (Discrete integration by parts formula). Let v,f,g be discrete functions such that 
vo = wjv+i = 0. Then we have the following identities: 

9{d+f) = - / (d^g)f + g N +ifN+i - <7o/o ; (2.17) 

[0,1) .7(0,1] 

9{dhf) = ( (mtg)(dif) - £flo(#/)o - ~gN + i{d^ f) N+ i ; (2.18) 

(0,1) "'[0,1) z z 

gv(d h v) = - [ \v\ 2 d h g + % f \dtv\ 2 d+g ; (2.19) 

(0,1) •'(0,1) Z -'[0,1) 



g{A h v) = - (d^v) {d£g) - (d£v) go + {d h v) N+1 g N+1 ; 

(0,1) J [0,1) 



(0,1) 2 J [0,1) 



Proof. Let us begin with (I2.17P : 

N , , , N JV 

/ gdtf = ^E*( : ^F A )=E^-E^ 

J [0,l) i=0 V ft / j=0 3=0 

JV+1 JV+1 

= 03-1/3 ~ £ + 3^+i/jv+i _ fo/o 

3=1 3=1 
JV+1 

= ~ h £ { 9l h ' 1 ) h + 9 N+1 f N+1 - 9ofo- 

3=1 

In order to prove (|2.18[) . using (J2TT2J , we do as follows: 

<A/ = \[f gdnf+l gdtf) 



(0,1) 



7 N , N 

3=1 3=1 

, JV-1 , JV 

3=0 3=1 

, JV h h 

2 + fl , J+i)0 9 fe/)j - 29a( d hf)o ~ -^9N+i{dif)i 

3=0 



(2.20) 



gv(A h v) = - f (d+v) 2 (m+g) + ± / M 2 A fe9 ; (2.21) 

(0,1) ^ [o,i) ^ -'(O.l) 

i r i 2 i i 2 

gA h vd h v = -- \d£v\ d£g+ - \(d^v) N +i\ ffjv+i - - I (5^ v)o I go. (2.22) 
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To prove (|2.19|) . using the fact that vq = hjv+i = 0, and successively ()2. 18[) . (|2.14p . (|2.15|) and 
(|2.17|l . we obtain: 



2 / gvd h v = 2 f m+(vg)(d+ 

J (0,1) "'[0,1) 



'[0,1) 



(m+v)(mtg) + ^(d+v)(d+g) ) (d+v) 



[ (m+g)d+(\v\ 2 ) + ^- f (d+v) 2 (d+g) 

J [0,1) 1 "'[0,1) 

(a-( m + ff ))H 2 + ^ f (d+vf(d+g) 

1,1] Z "'[0,1) 

{d h g)\v\ 2 + £ / (d+vf(d+g). 

1,1) Z "'[0,1) 



For (pT20]) . using pT7)l . we write 

/ g(A h v) = [ gd^d+v 

Jf0.11 JlO.l) 



(0,1) J (0,1) 

gdhdh v - £W+i((<9ftU)w+i - (9h «)jv+i) 

(0,1] 

= -/ u) (^h ff) - (d£v) g + (d^ v) N +igN+i- 

J [0,1) 

From (|2.20p . we prove (|2.2 1 1) . using wo = «w+i = and Lemma 12.61 
gvA h v = - f (d+v)(d+(gv)) 

(0,1) J [0,1) 



/ |9,>| 2 "i+g- / d+vm+vd+g 

J [0,1) "'[0,1) 
"'[0,1) Z "'[0,1) 



'[0,1) ^ -'(0,1) 



\g+\j \v\ 2 A h g. 

z -'(0,1) 

Finally, in order to prove (|2.22[) . we first remark that, using Lemma [2.61 
{A h v) {d h v) 3 = {d^(d+v)) 3 {m- h {div)) 3 = (|9+«| 2 ))j 
and therefore (|2.22|l follows from (|2.17|) : 

/ gA h vd h v = - [ gdl (|9^| 2 ) 

"'(0,1) z ■'(0,1) 

= \\ gdn {\dtv\ 2 )-lg N+ i{\(d+v) N+1 \ 2 -\(d~v) N+1 \ 2 ) 

Z J{0,1] z 

= ~h [ \dh v \ 2 dtg+ - v) N+1 \ 2 g N+1 - i \{d£v) \ 2 g . 
/ -'[o,i) £ 1 

This concludes the proof of Lemma 12.71 □ 
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2.4 Computation of the conjugate operator 

Set p = exp(— scp), ip given by (|2.1|) . and 

v(t,x) = p~ 1 (t,x)w(t,x) = e s,p{t ' x) w(t,x) and P h v := - (d tt - A h ) (pv). (2.23) 

P 

Proposition 2.8. The conjugate operator Ph can be expanded as follows: 

Phv = dttv + 2d t v h w 1 + -= Afc« — 2(4 -u v . (2.24) 

p p \ 2 p ) p p 

Proof. Identity (|2.24l) can be deduced easily by explicit computations based on Lemma 12.51 In- 
deed, 

P h v = - [dtt - Afc ](pv) = 9«w + 2<9tu h« — -. 

P P P P 

But, using (|2.16p . we get 

Ah(pv) * m h p 9 h p A h p 

= A h u h 2dj l t> h m h v . 

P P P P 

Besides, from (f2~TT|l . 

/ m h p\ /i 2 / A h p\ /i 2 
' ' = 1 + — and m h v = v + — A h u, 



p y 4 v p / 4 

which immediately yield identity (I2.24p . □ 

One step of the usual way to prove a Carleman estimate is to split the operator Ph into 
two operators Ph,i and Ph,i (detailed in Section [2.6[) , that, roughly speaking, corresponds to a 
decomposition into a self-adjoint part and a skew-adjoint one, and then to compute and estimate 
the scalar product 

/ P h ,iv P h , 2 v dt. 

-T J (0,1) 

But we first need to give a more precise expression of PhV, using the following equalities: 
Proposition 2.9. The coefficients in the expression of Ph can be expanded as follows: 

$t $i i 2 2 2 2 2 2 

= —sXipdtip, = s X tp {d t tp) — s\ <p (d t ip) — sXpdttip, (2.25) 

P P 

?hE = - sX A 1 , ^hE = s 2 x 2 A 2 -sX 2 A,-sXA i , (2.26) 
P P 

where the coefficients (Ai, A2, A3, A4) are given by 

g — sip(t,x + o-h) 

[pd.i>] (t, x + ah) e _ sip(t>x) da, (2.27) 
- \a\) [ V 2 (d x xb) 2 ] (t, x + ah f *^X*) h da > (2 - 28) 

t ,Xj -\-<rh) 

da, (2.29) 







Ax(t,x) = 








A 2 (t,x) = 


/> 










A 3 (t,x) = 












A 4 {t,x) = 





' e ~SLp{t,x) 

g—s(p(t.Xj-\-ah) 



e -.„(«,.) ^ ( 2 - 3 °) 
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Proof. Since p — e sv and ip = e x ^ , identities (I2.25|l are straightforward. 
Getting (|2.26[) is more technical. We write 

(d h p) j = Pj+1 2fe Pj ~ 1 = J ^ d x p{x) dx = ]^ j ^ d x p(x 3 + ah) da 

and since d x p = —sXp<pd x ip, we get (!2.26P i : 

^) (*) = — f IMmit + 'h) **'?*' 1 * da = -sXA hj (t)- 

Similarly, the proof of (|2.26[b relies on the usual Taylor formulas in integral form 

f(x ±h) = f(x) ± hf'{x) + f (1 =p ct)/"(x + ah) da. 

Jo 

Indeed, applying this identity to p, 

Pj±i = Pj ± hd x p(xj) + h 2 I (1 =p a)d xx p(xj + ah) da, 

Jo 

and therefore, 

(A h p) 3 = Pj+1 ~ 2 ff + ^ = ^ (1 - \a\)d xx p( Xj + ah) da. 

Since d xx p = s 2 X 2 pep 2 (d x tp) 2 — sX 2 p<p(d x ij}) 2 — sXptpd xx ip, we immediately deduce (|2.26[) c>, □ 

Remark 2.10. The coefficients of are intrinsically defined on the grid {j/i}je[i,JVl an d not 
for x £ [0, 1] as formulas ()2.27|) - (|2.30|) may imply. But it turns out that these formula induce 
a natural continuous extension of these coefficients that is easier to handle. We shall therefore 
identify these coefficients with their continuous extension given by (|2.27|) - (|2.30|) without confusion. 

2.5 Preliminary estimates 

Before going into the proof of the Carleman estimate itself, done in Section 12.61 we give here 
several key approximations on the coefficients Aj defined in (|2.27|) — (f2.28[) — (f2729|) — (|2.30[) and their 
derivatives. 

To begin with, we shall introduce the Landau notation 0\(e) to denote functions / = f(t, x) 
that satisfy, for some constant C independent of e > but that might depend on A, |/j < Ce. 
We are then in position to state the following basic estimates: 

Lemma 2.11. For all X > 0, s > and h > with sh < 1, for all a 6 [—2,2] and (t,x) £ 
[— T, T] x [0, 1] , 

p(t,x + «h) = e-°^ = 

p(t,x) e -s V (t, X ) \ I > V J 

a„f ""^, + r' l) Vo>(»'.); a.C ""'* + f ) 'Uo > (.t). (2.33) 



p(t,x) ) ' V p(P,x) 

Proof. Since the function ip is smooth and bounded on (— T, T) x (0, 1), we have ip(t, Xj + ah) — 
ip(t,Xj) + O(h) and therefore 

tp(t,x + ah) = e ^, X]) +xo(h) =e ^(*.^)( 1 + c , A (/ l )) = (p (t,x) + O x (h). 
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Therefore, we easily get (|2.3ip since 

p(t,x + ah) = e -"P<.**)+OxM = p (t,x)(l + Oi{sh)). 

Similarly, 

g / p(t, X + ah) \ _ q f-sip{t,x+ah) e sip(t,x)\ 

V P(t,x) J X \ ) 

= s\ [-(p(t, x + ah)d x ip(t, x + ah) + tp(t, x)d x ip(t, x)) ■ & _ sv{tiX) 
so that 

d x = sXO x (h)(l + Ox(sh)) = Ox(sh), 

which concludes the proof of (|2.32l) . left. 

Of course, other estimates in (|2.32[) - (|2.33[1 can be proved following the same ideas. Details 
are left to the reader. □ 

We can now give good approximations of the coefficients Aj: 
Lemma 2.12. Set 

fi = <pd x tp, fa = (d x ip) 2 , h = f(dxtp) 2 , f4 = tpdxxip- 

Using the notations Aj defined in Provosition \2.!A for (t,x) £ [— T,T] x [0, 1] and j £ {1,2,3,4}, 
we have: 

• On the 0th order derivation operators: 

Aj = fj + Ox(sh) = mt(Aj) + Ox{sh) = ml{Aj) + Ox{sh) = m h (Aj) + Ox(sh) ; (2.34) 

• On the 1st order derivation operators: 

d h Aj = d x f s + Ox{sh) = d+Aj + Ox(sh) = d^Aj + Ox{sh), (2.35) 
d t A j = dtf J +Ox(sh) ■ (2.36) 

• On the 2nd order derivation operators: 

AhAj = d xx fi + Ox{sh), dttAj = dufj + Ox(sh). (2.37) 
Proof. Let us first notice that that all the coefficients Aj can be written as 

pi g-s^fi.x+CTh) f 1/2 if j = 1 

Mt,*)= ]_HV)fi{t,x + oh) e _ tvM da, fe(a) = | (1 _ H) 

Using Lemma 12.111 and the regularity of %j) and <p, one can write 

- s<p(t,x-\-ah) 
-s<p(t,x) 

-1 



Aj{t,x) = j ^ fi, (a)fj (t, x + ah) £ ^ da 



N {a){f ] {t,x) + O x (h))(l + Ox(sh)) da = f 3 (t,x) + O x (sh). 

Let us remark that it also yields the same expansion for Aj(t, x + h) up to an error term of order 
Ox{h), from which one easily concludes (|2.34|l . 

For the first-order derivatives (|2,35[) . we can write 

1 f 1 

d h Aj(t,x) = - J dxAj(t,x + ah) da 

= 2 /_! W( ff ) a '^(*' a!+ (° + e -s V (t,x+ a h) d0ld<J 

/•! /•! /g— «lp(t,»+(a+o , }A) \ 

+ 2 J_ i J_ i ^ CT )fi^ x + { a + (T ) h ) a - [ e s V (t,x +ah) J dada - 
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But using 

1 f 1 f 1 / e ~s^{t,x+(a + a)h) \ 



Therefore, we only have to estimate 

^ r 1 r 1 e —s<p(t,x+(ct+<r)h) 

which can be done by using d x fj(t,x+ (a + a)h) = d x fj(t,x) + Ox(h), (|2.3ip and the fact that 
f^_ 1 fj,j(a)d<j=l. This yields 

(t, a;) = d x f j (t,x) + O x (sh). 

Of course, similar computations can be done for d^Aj, Aj. 

For (|2.36p . the idea is the same: we use the integral expression of the coefficients, check that 
the derivatives from the ratio of exponentials are of order 0\(sh) and can therefore be neglected 
due to (|2,32[) . and then proceed as above. Details are left to the reader. 

For the estimates on the second order derivatives (|2.37p , this proof applies again and is there- 
fore omitted, using this time the second order estimates (|2.33p . □ 

To summarize the results detailed in Lemma [2.121 we have proved that 

Ai ~ (pd x t/>, A 2 ~ tp 2 (d x ip) 2 , A 3 ~ ip(d x i()) 2 , A& ~ <pd xx ip, 

up to error terms in (D\(sh), and these expressions can be differentiate twice, still with an error 
term of the order of Ox(sh). 



A more precise expression of PhV can now be deduced from Propositions 12.8 1 and 12.91 

Phv = dttv — 2s\ifidtipd t v + s 2 \ 2 ip 2 (d t ip) 2 v — s\ 2 (p (dttp) 2 v — s\ip{dttip)v 

- ^1 + y (s 2 X 2 A 2 - sX 2 A 3 - sXA 4 )^ A h v + 2sXA 1 d h v ~ {s 2 X 2 A 2 ~ sX 2 A 3 - sXA 4 )v. 

In order to simplify the notations, we set 

h 2 

A = — (s 2 X 2 A 2 - sX 2 A 3 - sXA 4 ), 

so that Ph can be rewritten as 

Phv = dttv — 2sXtpdt^dtv + s 2 X 2 (f 2 (dt*l>) 2 v — sX 2 tp (d t *l>) 2 v — sXtp(dttip)v 
- (1 + Ao) A h v + 2sXA 1 d h v - (s 2 X 2 A 2 - sX 2 A 3 - sXA 4 )v. 

Note that Ao is expected to be small. Indeed, using Lemma 12.121 one easily gets that Ao is 
in (D\(sh) and that the same holds true for the following expressions: 

Ao,rn h A , m^Ao, d h A , d^A , d t A , A h A ,dttA all are G\(sh). (2.38) 

We emphasize that this term Ao is a purely numerical artifact which does not have any continuous 
counterpart. 
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2.6 Proof of the Carleman estimate 

In this section, we focus on the proof of the discrete Carleman estimate (|2.7|) given in Theorem l2.3l 
We first set 

Ph.iv = duv - A h v(l + A ) + s 2 X 2 [cp 2 {dti>f - A 2 ] v , (2.39) 

Ph,2v = -sX 2 [<p\d t ^\ 2 -A 3 ]v- 2s\ [ipd t ipdtv - A 1 d h v] , (2.40) 

R h v = s\ [<pd tt ip - A A ] v , (2.41) 

so that we have Ph,iv + Ph,2V = PhV + RhV, and in particular, 

r-T r r T r 

\P h v + R h v\ 2 dt = / / \P hA v\ 2 dt + 



h,2 l 

tJ(o.i) J -t J (0,1) J-tJ(o,i) 

+ 2 f f P hA vP h , 2 vdt. (2.42) 

J-TJ(0,1) 

We will then prove the following: 

Lemma 2.13. There exist A > 0, so > 0, eo > and a constant Mo > such that for all 
s G (so,eo/h), for all v satisfying vo = wjv+i = and v(±T) = dtv(±T) — 0, 

T f \d t v\ 2 dt + s r / \d+v\ 2 dt+ s 3 r / \v\ 2 dt 

-T J (0,1) J-T J[0,1) J -T J (0,1) 

+ [ [ \P h ,!v\ 2 dt+ [ [ \P h , 2 v\ 2 dt < M [ [ \P h v\ 2 dt (2.43) 

J-T J (0,1) J -T J (0,1) J-T J (0,1) 

+M s f \(d^v) N+1 \ 2 dt + Mos [ [ \hdid t v\ 2 dt. 

J-T J-tJ[0,1) 

Proof. We will begin with calculating and bounding from below the product 

rT 



Ph,lV Ph,2V dt. 

T J (0,1) 

Step 1. Explicit computations of the cross product. 

The proof of estimate (|2.43|) relies first of all on the computation of the multiplication of each 
term of Ph,iv by each term of Ph,2V. We write 



fT r 3 2 

/ / Ph,lvP h ,2Vdt = J2J2 1 ^ 
J-T J (0,1) i= l i=1 



where Uj denotes the product between the i-th term of Ph.i in (12.391) and the jr-th term of Pk,2 
in (|2.40[) . We now perform the computation of each ly term. 

Of course, we shall strongly use below the properties of v on the boundary and in particular 
that «(±T) = dtv(±T) = 0, Vo(t) = VN+i(t) — and also the fact that dtVo(t) = dtVN+i(t) = 
for all t € (— T, T). 

We shall also use the results of Lemma 12.121 which will allow us to simplify most of the 
expression in which the coefficients Aj appear. We recall to the reader that we will use the 
notation (|2.3p for the discrete integrals in space. In order to simplify notations, we will also omit 
u dt" in the integrals in time (that are continuous ones). Therefore: 

f T I f = hf^[ T f 3 {t)dt and f T f f = hjr f T fd(t)dt. 

J-T J (0,1) j=1 J-T J-T J [0,1) j=0 J-T 
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Computation of In. Integrating by parts in time, 

In = - sA 2 / / d t tv((p\d t ip\ 2 - A 3 )v 

J-T J (0,1) 



J-T J (0,1) 1 J -T J (0,1) 



-TJ(0,1) * J -T J (0,1) 

S A 2 r f \d t v\ 2 <p(\dti>\ 2 - i&^i 2 ) 



sX 
2 



TJ(0,1) 
2 /*T 

v| 2 C>tt(tp|9t'!/'| 2 — tpl^VI 2 ) 

T J (0,1) 



+ s/ / Ox(sh)\d t v\ 2 + S [ [ Ox(sh)\v\' 

J-T J (0,1) J-T J (0,1) 

using A 3 = ^l&VI 2 + A (s/i) and d tt A 3 = 9„ (^Vj 2 ) + Ca (sft). 
Computation of In ■ 



hi = - 2sA 



/ / duv(<pdtipdtv — Aidhv) 

J-T J (0,1) 



fT r 

s\ I / \d t v\ 2 ipd t ti> + s\ 2 

-T J (0,1) 



2sA 



T 

-T J (0,1) 



But, by ([2~19j) 

fT 

2s\ 



<-TJ(0,l) J -T J (0,1) 

Therefore, using Lemma 12.121 for dhAi and dtAi, we get 

fT 

1 12 = s\ 

' -T J (0,1) J-T J (0,1) 

rT r „\ rT 



f 7 


\d t v \ 2 ip \d t ip 2 


./-T ^(0,1) 








2sA /" 


/ Aidtvdhdtv 


J-T J 


(0,1) 


2 d h Ai - 






z J-TJ(0,1) 



/ / \d t v\ 2 ip{d tt ^ + d xx ^) + s\ 2 f f \d t v\ 2 v{\dtM 2 + \9^\ 2 ) 

J-T J (0,1) J —T J (0,1) 

2s\ 2 [ T [ <p{d t v){d h v) d t ip d x ip- ^ / T / ifea+ft^^+Ax 

J-T J (0,1) * J-TJ\0,1) 



T J (0,1) z J-TJ[0,1) 

O x (sh)\d t v\ 2 + s f T f O x {sh)d t vd h v. 

T J (0,1) J -T J (0,1) 



Computation of hi- Using ([2~^T]| and l|2.38p . 

J-T J (0,1) 



T / l^l'm+CCl + Ao)^!^! 2 -^)) 

-T J [0,1) 



2 /-T 



+ ^-/ / |<A fe ((i + A,)(¥Wf-A 3 )) 

^ J -T J {0,1) 

f T f \d+v\ 2 v (\dti>\ 2 -\d x ip\ 2 ) 

J-T J[0,1) 
^ J-TJ(0,1) 

+ 8 f I Ox(sh)\d+v\ 2 + S f f Ox(sh)\v\ 
J-T J [0,1) J —T J (0,1) 
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We do not develop the term in da; X {ip\dtip\ 2 — ip\d x tp\ 2 ) since it is uniformly bounded with respect 
to s. 

Computation of I22 ■ We can split this term in two parts as follows 



I22 = 2sX 



/ f A h v(l + A )tpd t ipdtv~2s\ [ [ A h v(l + A )A 1 d h v. 
J-tJ(o,X) J -t J (0,1) 



l22a ^22b 

To compute /22a, we use A;, = d^d£, dtvo = dtVN+i = 0, and formula (12.15[l and (|2.17|) : 



/22a = -2s\ [ f (d£v)d£(0-+Ao)<pdtil)dtv) 

J-T J [0,1) 

= -2s\[ T [ (d+v)m+(d t v)d+((l + A )ipd t i>) 

J-T J[0,1) 
rT 



+ sx[ f \d+v\ 2 d t mt{{l + A )^). 

J-T J[0,1) 

-2s\ [ T [ (d+v)m+(d t v)d x (<pd t i>) + s\ f T f \d+v\ 2 d t Wdt^) 

J-T J [0,1) J-T J [0,1) 

+ sf T [ O x (sh)\d+v\ 2 + s f T f O x (sh)(d+v)m+(d t v). 
J — T J [0,1) J-tJ[o,d 



-TJ[0,1) J —T J [0,1) 

For the computation of /22b, we rather use (|2.22[) and (12,381) 

rT 



122b = — sX 



I I \d+v\ 2 d+{{l+Ao)A 1 ) 

J-T J [0,1) 

+ s\ j T ((1 + A )A 1 )(t, l)\(d^v) N+1 \ 2 

rT 



sX [ ((l + A )A 1 )(t,0)\(d+v) \ 2 

J-T 

sA [ T [ \d+v\ 2 d x {ipd x i>) + s f T f Ox(sh)\d+v" 2 

J-t J ro.ii J-T J \0.1) 



TJ[0,1) J-T J [0,1) 

+ .sA / T ([^](t, 1) + Ox{sh)) \(d^v) N+1 \ 2 

-T 

•sA f T (&d x ^](t,0) + Ox(sh))\(d+v) \ 2 . 

-T 
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Therefore, 7 22 = 7 22a - hib gives 



I22 = — 2sA 2 / / (dhv)m^(dtv)tpd x ipdtip 
J-T J [0,1) 

f T l \d£vW{\d*i>\* + |a,vi 2 ) 

7-T J[0,1) 



+ sA 2 



T J [0,1) 

T 

+ ».|2, 



+ sX / ip(d u ip + d xx ip) 

J-T J[o,i) 

- s a y ([^.vK*, 1) + o x {sh)) \(d^v) N +i\ 2 

+ s\ ({<pd x i>}(t,0) + Ox(sh)) \(d£v)o\ 2 

+ s[ f Ox(sh)\d+v\ 2 + s [ [ Ox(sh)(d+v)m+(d t v). 

J-TJ[0,T) J-T J [0,1) 

Computation 0//31. From Lemma [2, 121 one easily obtains: 

7 3 i= -s 3 X A [ [ \v\ 2 (ip 2 \d t ij\ 2 - A 2 )(<p\dtil>\ 2 - A 3 ) 

J-T J (0,1) 

= - s 3 a 4 / f |»iV 3 (ia^i 2 -iwf) ! + « 3 ( T I o x {sh)\v\ 2 . 

J-T J (0,1) J —T J (0,1) 

Computation 0/732. Finally, from (|2.19|l of Lemma 12.61 we get 

rT 



7 32 = - 2s 3 A 3 / / v{ip 2 \d t ip\ 2 - A 2 )(<pd t ipd t v - Aid h v) 

J-T J (0,1) 

= s 3 A 3 f T f \v\ 2 d t ({<p 2 \dti}\ 2 -A 2 )<pd t ij) 

J-T J (0,1) 



-T./(0,1) 

-s 3 X 3 f T [ \v\ 2 d h {Ai(<p 2 \d t ip\ 2 - A 2 )) 

J-T J (0,1) 

+ T f \hdtv\ 2 dt{Ai{<p 2 \dtM 2 - M)). 



T J (0,1) 



But, according to Lemma 12.121 

dt{^ 2 \dti>\ 2 - A 2 )<pdti>) = 3\p 3 \dt*p\ 2 (l&VI 2 - \dxi>\ 2 ) 

+ Z^dtt^dt^ 2 - ip 3 d u il)\d x il>\ 2 + Ox(sh) 

and 

dhiA^cp 2 ^] 2 - A 2 )) = 3\cp 3 \d x 4>\ 2 (\dti>\ 2 - |a^l 2 ) 

+ <p 3 d xx tp\d t ^>\ 2 - S(p \d x ip\ 2 d xx ip + 0\{sh). 

Thus we obtain 



T 

T J (0,1) 



7 32 = 3s 3 a 4 / / MvaW-ia^l" 



+ 3s 3 A 3 / / |d| 2 (^ 3 (|9*V'| 2 9ttV' + I^VI 2 ^^^) 

i-T i(0,l) 

— S 3 A 3 / / |^| 2 <(5 3 (|c?a;'0| 2 C?tt'0 + |(9tV , | 2 ft I;:E l/>) 

7-T i(0,l) 

+ s/ / Ox(sh)\d+v\ 2 + s s f f Ox(sh)\v[ 

J —T J (0,1) J-T J (0,1) 
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Final computation. Gathering all the terms, one can write 

Ph,lVP h ,2V = I v + Id t ,Vv + I{0,1} + ^Tych, (2.44) 



J-T J(0, 



,1) 

where I v contains all the terms in \v\ 2 : 



T 



Iv = I I H 2 F, (2.45) 

T J (0,1) 



with F given by 



\2 „\2 



sX 

F = 2~dtt(<p\dtip\ 2 - ^I^VI^) + ^-d xx (ip\d t tk\' 1 - ^S^l 

+ s 3 X 3 ip 3 (\d t ip\ 2 - \d x i)\ 2 ){d t ti> - d xx i>) + 2s 3 X 3 <p 3 (\d t ip\ 2 dttip + \d x i>\ 2 d xx il>) 
+ 2s 3 \ 4 <p 3 (\d t ip 2 - \d x ip\ 2 ) 2 + s 3 O x (sh) 
= s 3 X 3 ip 3 (\dtip\ 2 - \d x t[>\ 2 )(dttip - dxxip) + 2s 3 \ 3 ip 3 (\dtip\ 2 dttip + \d x t[>\ 2 d xx ip) 
+ 2s 3 \ 4 p 3 (\d t iP\ 2 - \d x iP\ 2 ) 2 + s 3 O x (sh) + sO x (l) ; 
Id t ,Vv contains all the terms involving first order derivatives of v. 

I df . Vv =2s\ 2 ^ f l»-.l 2 ..l».'.l 2 ■ o.x2 F f ,2 



/ \d t v\ 2 <p\d t ip\ 2 + 2s\ 2 [ [ \d+v\ 2 <p\8 x ip\ 2 

J (0,1) J-T J [0,1) 

-2sX 2 / ip{d t v){d h v) dti> d x tp-2sX 2 / / ip{&Z 

J-T J (0,1) J-T J [0,1) 

+ sX / \dtv\ 2 ip{dttip + dxx^p) + s\ I I v\ 2 (p(d tt ip + d xx tp) 

J-T J (0,1) J-T J [0,1) 

-A 

J-T J (0,1) 

+ S / T / (Oa(^)|^«| 2 + OaM™* 

J-T -/[O,!) V 7 

/{o,i} contains all the boundary terms: 



! / T / (oA(^)|a^| 2 + OA(^)a t ^ h «) 

J-T J (0,1) V 7 



/{o,i } = - s\J ([<pd x i/>](t, i) + Ox(sh)) |(a ; ;«)iv +1 j 2 

+ S A y {[<pd x iP](t,0) + Ox(sh))\(d+v) \ 2 ; 
^Tych contains the term corresponding to the Tychonoff regularization: 

lTy* = ~[ T [ Ihd+dtvfd+A,. (2.46) 
z J-tJ[o,i) 

Step 2. Bounding each term from below. 

In the sequel, M > and C > will denote generic constants that all are independent of h and s 
but may depend on A. 

Step 2.1. Dealing with the order terms in v. 
We have: 



F = s 3 X 3 <p 3 G + s 3 O x (sh) + aO x (l), 



with 



G = 2\(\d t ^\ 2 - \0xip\ 2 ) 2 + (\dtip\ 2 - \dxip\ 2 )(d tt *l> - d xx ii>) + 2(\d t *p\ 2 dui> + \dxip\ 2 dxx*P) 

= 2\(\d t ib\ 2 - \d x *p\ 2 ) 2 + (IdtVI 2 - \d x ip\ 2 )(d tt *p - dxx*P) + 2(\d t ib\ 2 - \dxM 2 )dtt*P 
+ 2\d x ip\ 2 (d„tp + d xx ip). 
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But dutp + 9 xx ip = 2(1 — P) > and inf( ,i) I^VI 2 = 4inf( 01 ) \x — x a \ 2 is strictly positive since 
x° i [0, 1]. Therefore, setting X = \dtiti\ 2 - \d x ^\ 2 , we have 

G > 2 AX 2 - 2X(3/3 + 1) + c, with c = 16(1 - /3) inf |x - x°| 2 > 0. 

Thus, there exists A > large enough such that for all A > A , this expression can be made 
strictly positive. Therefore, for A > Ao, we get a positive constant c* > independent of A such 
that 

I v >2c,s i \ 3 f f cp 3 \v\ 2 -(s 3 Ox(sh) + sOx(l)) f f \v\ 2 . 

J-T J (0,1) J-T J (0,1) 

Thereby, bounding ip from below by 1, we can choose so (A) such that for all s > so (A), 



I v > c*s 



3 A 3 f T f \v\ 2 -s :i O,(sh) f 1 f \v\ 2 . (2.47) 

J -T J (0,1) J-T J (0,1) 



From then on, we fix A > Ao. 

Step 2.2. Dealing with the first-order derivatives. 

The idea is to show that the terms in which sA 2 appears are positive up to a small error term, 
and then check that the terms in sA are strictly positive. 
On the one hand, 



/ / <p(dtv)(dhv) <9tV < \ I I v\dtv\ 2 \dt^\ 2 + \ [ [ v\dhv\ 2 \d x i>\ 2 

J-T J (0,1) ^ J-T J (0,1) i J-T J (0,1) 

< \ r [ <p\d t v\ 2 \d t ip\ 2 +t- r f ^\d+v\ 2 \d X i>\ 2 +o X (sh) r / \d+v\ 2 

z J-T J (0,1) A J-TJ\0.1) J-TJI0.1) 



and similarly, 



/ / ^(3h«) m+(d t v) dttp d x tp 

J-T J{0,1) 



< If I ^\dtv\ 2 \d x ^\ 2 + \ f f V \m + h (d t v)\ 2 \d t i>\ 2 

* J-T J[0,1) 1 J-T J\0,1) 

< \ r f v\d+ v \ 2 \d^\ 2 +i r [ <p\d t v\ 2 \d t i>\ 2 +o X (sh) r / \d tv \ 2 . 

z J-tJ[o,i) z J-tJ(o,i) J -t J (o,i) 

Therefore, 

2sA 2 f f \d t v\ 2 <p\d t ^\ 2 + 2s\ 2 f I \d+v\ 2 ^\d x ^\ 2 

J-T J (0,1) J-T J [0,1) 

-2sA 2 / / <p(d t v){d h v) dtip d x ip - 2s\ 2 / ip{dZ 

J-T J (0,1) J-T J [0,1) 

>-so X (sh) r f \d+v\ 2 -so X (sh) r [ \a, 

j-tJiod J -T J (0,1) 



<tv\ 2 . 



On the other hand, focusing on the terms in sA, we have dtttp + d xx ip = 2(1 — /3) > 0, and 
then, bounding ip — e A * from below by 1, we obtain: 

7 9t ,v, > 2sA(l - 0) [ f \d+v\ 2 +2s\(l- /3) [ [ \8 t v\ 2 

J-T J [0,1) J-T J (0,1) 

-sO x (sh)f T f \d+v\ 2 - sOx(sh) f T [ \d t v\ 2 , (2.48) 

J-T J [0,1) J-T J (0,1) 
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where we used that, by Cauchy Schwartz, 

T / ^ < r [ \8+v\ 2 , r f \m+(d t v)\ 2 < r f la^. 

■T J (0,1) J -T J [0,1) J -T J [0,1) J —T J (0,1) 

Step 2.3. The boundary terms. 

Since min.(_T,T)x(o,i) vd^ip > (recall x° < 0), then there exists ei(A) > such that taking 
sh < £i(A), 

0\(sh) < min lv(i, x)d x ib(t, x)\ . 

(t,x)e(-T,T)x(0,l) 

Therefore, 



{ o,i } =-sA / ([^](t,i) + o>(sh))\(d l ;v) N+1 \ 2 + s \ f (&d X Tp](t,o) + o4sh))\(d+v) \ 2 

J — T J-T 

>-2sX [ipd x i>](t,l)\{d^v) N+1 \ 2 >-sCx |(<9,»iv + i| 2 . (2.49) 

J-T J —T 



Step 2.4- The Tychonoff regularization. 

Let us recall that d^Ai = \ip\d x i>\ 2 + tpd xx ip + Ox(sh) = O x (l) since sh < e(A). Thus 



hych = ~ f l \hd+d t v\ 2 d+A 1 > -sCx [ [ \hd+d t v\ 2 . (2.50) 

2 J-TJ[0,1) J-tJ[0,1) 

|2 



Noticing that \(p\d x i/}\ + ifid xx i/j > 0, /Tych < and cannot be made positive. This is not only a 
technical matter, since otherwise we would get uniform observability results for the semidiscrete 
wave equations, which is proved not to hold in [25] . 

Step 3. Proof of Lemma [27131 

Recall that A is fixed from Step 2.1. 

Collecting the results (|2.48|l - (|2.50|l . of Step 2, we have proved that for s > so(A) and sh < £i(A), 
/ P h ,ivP h , 2 v > 2sX(l-/3)F f \0+v\ 2 + 2sA(l-/3) F f \d t v\ 2 

T J (0,1) J -T J [0,1) J —T J (0,1) 

-sOx(sh) [ f \d+v\ 2 - sOx(sh) f I \d t v\ 2 
J-tJ[o,i) J-t J (0,1) 

+c„ S 3 A 3 F f \v\ 2 - s 3 Ox(sh) F f \v\ 2 

J-T J (0,1) J-T J (0,1) 

-sC x F \(d^v) N+1 \ 2 - sC x F f \hd+0 t v\ 2 . 

J-T J-T J [0,1) 



Therefore, taking sh small enough such that 

f e A 3 

Ox(sh) <min|A(l-/3) J ^-,ei(A) 

which defines e(A) > 0, we obtain, for some constant Mi — Mi (A) > 0, 
\d t v\ 2 + sF f \d+v\ 2 + s 3 F 

T J (0,1) J-T J [0,1) J -T J (0,1) 



r T f rT r T 

<Mi / Pi,hvP 2 ,hV + Mis 

J -T J (0,1) J-T J-TJ[0,X) 



if \{d^v) N+1 \ 2 + M lS f f \hd+d t v\ 2 . (2.51) 

J-T J-T J 10.1) 
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Now, we use that from (I2.42|l , 

2 I f Pi, h vP2, h v+ f I \Pi.hv\ 2 + I I \Pi,hv\ 2 



T J (0,1) J -T J (0,1) J -T J (0,1) 

■T r fT 



<2 f f \P h v\ 2 + 2f f \R h v\ 2 
J-tJ(o,i) J -T J (0,1) 



1(0,1) J -T J (0,1) 

where R^v is given by (12.41[) . which yields, for some M2(A) > 0, 
T I \d t v\ 2 + sf T f \d+v\ 2 + s 3 f T f \v\ 2 

T J (0,1) V -TV [0,1) J — T J (0,1) 

+ [ f \Ph,iv\ 2 + f f \Ph,2v\ 2 < M 2 [ [ \P h v\ 2 



-T J (0,1) J -T J (0,1) J -T J (0,1) 



r l \R h v\ 2 +M 2 s r \{d^v) N+1 \ 2 +M 2 s r f i^t^i 2 . (2.52) 

J -T J (0,1) J-T J-T J [0,1) 



+ M: 
Therefore, since 

r f \R h v\ 2 <s 2 x 2 r f \v\ 2 {< P d tt 4>-A i ) 2 < s 2 o X {i)f \v\ 2 , 

J -T J (0,1) J-T J (0,1) J (0,1) 

this term can be absorbed by the left hand side of (|2,52[l by taking s large enough, thus yielding 
to J233J). □ 

Proof of Theorem \2.3i The Carleman estimate (|2.7[) of Theorem 12.31 will now be deduced from 
Lemma l2.13l Actually, it simply consists in writing (12.43|) in terms of w instead of v, using that 
w — uexp(— sip) and, by construction, see (12.23[) . exp(scp)Lhiv — Ph,v. 
In particular, we have 

e^dtw = e aif 'dt{ve~ a,fi ) = d t v — sXvipdtii, 

e^d+w = e'ra+ive-**) = d+v(e sv m+{e-^)) ~ e*d+{e-^)m+v, 
and, since, similarly as in Lemma l2.11l 

e av m+{e- sv ) = 1 + O x (sh) and e s,fi d+ (e~ s ^) = -s\tpd x tjj + O x (sh), 
we get, for sh small enough, 

e 2sip \d t w\ 2 < 2\d t v\ 2 +2s 2 \ 2 ip 2 \d t 'ip\ 2 \v\ 2 , 
e 2s *\d+w\ 2 < 3\d+vf + 3 S 2 AV|d,V|>+«| 2 . 
Therefore, there exists Ms(A) such that 

e 2sv \d t w\ 2 +s f T f e 2sv \dtw\ 2 +s 3 C f e 2sv \w\ 2 

) J-tJ[0,1) J-T J (0,1) 

<M s s[ [ \d t v\ 2 + M 3 s[ [ \d£v\ 2 +M 3 s 3 f [ \v\ 2 . (2.53) 

J-T J (0,1) J-T J \0,1) J-T J (0,1) 



Now, it remains to estimate the right hand side of (|2.43l) in terms of w. For the boundary 
term, we use the fact that <p(t, •) is increasing and therefore 



| v N 
I h 



<\WN\ M t,l) = l{d - w)N+1 \ e ^,l)_ (2 54) 
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Finally, we bound the term corresponding to the Tychonoff regularization, using (|2.15p of Lemma| 
and d t v = e sv d t w + wdt(e 3lp ): 

hd+d t v = hd+(d t w)m+(e a n + h(d+w)m+d t (e a n 

+ rn+(d t w)hd+(e s n + (m+w)hd+ (d t e**). 

Again, similarly to Lemma l2.11l one can prove 

^ ( a t( e^)) = !*ten =0xiah)t and 

and deduce a bound of the form: 

\hd+8 t v\ < \hd£{d t w)\e sv {l + O x {sh)) + Ox(sh)\d+w\e sip 
+ O x (sh)\m+(d t w)\e s * + sO x (sh)\m+w\e s *. 

Therefore, one gets: 



.s 



hd+d t v\ 2 

T J [0,1) 

T c ,-T 



<M 4 sf f e 2sv \hd+d t w\ 2 +M 4 sO x (sh) f f e 2sv (\8+w\ 2 + \m+(d t w)\ 2 
J-tJ[o,d J-tJ[o,d 



-TJ[0,1) J-tJ[o,i) 
+ M 4 s s Ox(sh) f T f e 2ev \m+w\ 2 

J -T J[0,l) 



/•T 

< M 5 s 



f f e 2s *\hd+d t w\ 2 +M 5 aO x {8h) [ f e 2s *\d+w\ 2 (2.55) 

J-T J [0,1) J-T J [0,1) 

+ M 5 sOx(sh) f T f e 2s *\d t w\ 2 + M 5 s 3 Ox(sh) C f e 2sip \w\ 2 . 

J-T J (0,1) J-T J (0,1) 

Hence, combining (12.53p - (|2,54p - (|2.55|l . plugging them in (|2.43l) and choosing ah small enough 
so that MoMs,M^Ox{sh) < 1/2 and sh < so (given by Lemma |2.13[) . we obtain the desired 
Carleman estimate (|2.7p . □ 

3 Uniform stability estimates 

In this section, we state and prove uniform stability results for the semi-discrete framework, 
announced by (II. lip in the introduction. 

3.1 Statements of the results 

Similarly to Theorem ll.il we will prove the following local stability result: 

Theorem 3.1. Let m > 0, K > 0, r > 0, T > 1 and p h G L^< m (0, 1). 
Consider the equation 

( dtty 3 ,h - (A h y h ). + Pi,hVj,h = 9j,h, t G (0, T), j G [1, NJ, 
{ Vo,h(t) = Bh(t), VN+iAt) = gl(t), te(o,r), 
I VjAO) = yl h , d tVj , h (0) = vl h , j G [1, N], 

and assume that 

II?/MIIh1(0,T;L-(0,1)) < K 

and 

^W yU - r - (3 - 2) 
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Then there exists a constant C = C(T,m, K,r) > independent of h such that for all G 
Z/£° <m (0, 1), the uniform stability estimate holds: 

hh-Ph\\ L 2 i0tl) < C\\dt(d^y h ) N+1 [q h ]-dt(d^y h )N+i]ph]\\ L 2 {0iT) (3.3) 
+ C \\hd~l duyhWh] - hd hdttyh[Ph]\\ L2{0)T . L 2 h[0>1)) ■ 



Before giving the proof of Theorem [371] at the end of this section, we will begin by a stability 
theorem for the following inverse source problem, using the discrete Carleman estimate obtained 
in the previous section: 

Theorem 3.2. Let m > 0, K > 0, r > 0, T > 1. 

Let f h G Ll(0, 1), R h G H 1 (0,T;L^°(0, 1)) such that 

\\ R h\\ hHo,T;L^ (o,i)) < K and inf \R jA (0)\>r. (3.4) 
Let qn G Z/^< m (0, 1) and consider the semi-discrete wave equation 

( duUj, h - (A h u h ) j + q jth u jth = fj, h Rj, h (t), t G (0, T),j G [1, NJ, 

{ u o , h (t) = 0, u N+hh (t) = 0, t£(0,T), (3.5) 

[ u jth (0) = 0, d tUj , h (0)=0, j 6 jl, NJ. 

Then there exists a constant C = C(T, m, K, r) > independent of h and such that 

\\dt(dhU h ) N +i\\ L2{0T) + \\ h dt,dttu h \\ L2{0 T . L z {[01))) < C||/h|| L 2 (0jl) , (3.6) 

ll/fc|lz,»(0,l) - C \\ dt ( d h U h)N+\\\ L2(0 T) + C \\ hd h d ttUh\\ L 2 (0}T . L 2 h{[0Am ■ (3-7) 

Theorem 13.11 will then be a simple consequence of Theorem 13.21 see its proof in Section T3. 31 



3.2 Stability for the inverse source problem 



Before going into the proof of Theorem l3.2l we first recall some counterparts of the classical energy 
estimates for the solutions of the continuous wave equation in the context of the semi-discrete 
wave equation: 

Lemma 3.3. Let m > and q^ G Z/£°< m (0, 1). 

Let g h G L 1 (0, T; L\(0, 1)) and ( zl, z^) be discrete functions and let Zh be the solution of 

( d t tZ jth -(A h z h ) j +q j , h Zj,h=9jAt), * G (0,T),j G [l,iV], 

I z 0ih (t) = o, z N+1 . h (t) = o, te(o,T), 

{ Zj , h (0) = zl H , d t z jth (0) = zl h , j G [1, N\. 

Then, setting 

E'm = ||^"«fc(*)||l» ([0 ,i )) + Hft*h(t)llia (0 ,i) + lkfc(*)ll^(o,i) ' ( 3 ' 8 ) 
there exists a constant C = C(T,m) > independent of h and such for all t G (0, T), 

Eh(f) < C (j\dh Z °\\ L 2 aoA)) + |K|| L 2 ((M) + Wlz,l(0,T;££(0,l))) • ( 3 - 9 ) 
We also have the following "hidden regularity" property: 

\\(9hZ h ) N+1 \\ 2 L2{0 T) < C (Ph Zh\\ 2 L 2 h([01)) + \\zl\\ 2 L 2 i0A) + ll<?h|l!l(0,T;,L2(0,l))) ■ ( 3 -10) 
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Proof. The proof is somewhat classical, except perhaps for 1)3,101) . The quantity E^ is usually 
called the discrete energy of the solution Zh- We sketch it for the convenience of the reader since 
it will be useful in the sequel. 

Differentiating with respect to the time t, we obtain 



dE, 

~~dt 



z f 

-(i) = 2 / (dttZh — A h z h + z h ) d t Zh 

Ao,i) 

/ \g h (t)d t z h \+2{m + l) [ \z h d t z h \ 
J(o,i) ■'(0,1) 

< 2 ( Jg h (t)\ 2 ^ jE*ffl + (m+l)EZ(t). 



Therefore, 



d ^fM f , ,a\ (m + 1) — 

x < ( / \9h(t)\ 2 ) + 1 T VM 

'(0,1) 



dt 

and GronwalPs estimate then yields a constant C(T,m) such that for all t £ (0, T), 

El(t) < C(E* h (0) + IIShllino.T^o.u)), (3-11) 
which implies (|3.9[) providing a discrete Poincare estimate proved hereafter: 

N I j-1 

\z h \ 2 = hJ2 hJ2d£(\z h \ 2 ) k 
(o- 1 ) j=i \ k =o 



[ \dt 
J\o,i) 



\ 1/2 / \ 1, 

\dh z h\ 2 \ [] Qi \ m h z ^ 2 

1/2 / x 1/2 

I" 1 < / V(o,i) 



<2| / fel 2 / l*fc| a 



which implies 



M 2 <4/ fe| 2 . (3.12) 

(0,1) i[0,l) 

Therefore (pTT) implies (pi))) . 

Finally, to prove (|3.10[) . we use a multiplier type argument. Multiplying the equation of zh 
by j( z j+i,h — Zj-i,h) (which is a discrete version of xd x z), summing in j and integrating in time, 
we get (cf [251 Lemma 2.2] or the proof of (|4.20l) given hereafter in a more intricate case): 

N r T f T f \t 

/ / \d+ Zh \ 2 + x h (t)\ 

j=0 J ° 



J [0,1) 



2 I / g h xd h z h = J \(d h z h ) N+1 \ 2 (3.13) 
o J(o,i) Jo 



C h (t) = 2 xd h z h (t)d t z h (t). 
i(o,i) 



where 

-v,,cn - 

'(0,1) 

Of course, since each term in (|3.13l) is easily bounded by supp T j El(t) except for the term 
involving which can be bounded by ||Sh|| i i( i 2) sup[ T j yEfJt), we immediately obtain (|3.10|) 
from El). □ 
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Proof of Theorem \3.'A Step 1. Energy estimates. Set z% = d t Uh- Then, using the notation 
Lh[qh] = dtt — A h + q h , z h satisfies 

( L h [q h ]z h = f h d t R h , te (0,T), 

I z 0th {t) = z N+lih (t) = 0, te(0,T), (3.14) 

{ z h (0)=0, d t z h (0) = f h R h (0). 

We can apply Lemma 13.31 to Zh solution of (13.14[) since d t R h belongs to ^(0^-^1(0, 1)) and 
fh G L h (0, 1). In particular, if Eh denotes the energy of Zh (see (|3.8p l. we obtain, for all t G (0, T), 

Eh{t) < C {\\fhRh{Q)\\ 2 L 2 h(01) + ||/ft9t-Rft|lil( ,T;£2(0,l))) 

< C Il/h|li2(0,l) (\\Rh(0)\\ 2 L ^ i0tl) + \\Rh\\ 2 H l(0,T;L<x>(0,l))) ^ GK * Wfa II (0,1) ' ( 3 ' 15 ) 

where we have used that fl^O, T; L^(0, 1)) embeds into C([0, T]; L£°(0, 1)). 
Moreover, Lemma 13.31 also yields (|3.6p . Indeed, estimate p.lOp becomes here 

\\{d h Z h )jv + ij| L2(f)/r) + ||<9tZfe||i<*>(0,T;Z,2(0,l)) - C*-^ 2 11/^11^2(0,1) ' 

but Zh = 9t?ih and the operator hd£ is bounded uniformly in h. 

Step 2. The choice of the Carleman weight. Since we assumed T > 1, there exists x° < 
such that 

T > sup | as — x°\ [ = 1 + |a 

ze(o,i) v 

Therefore, we can choose /3 G (0, 1) and rj > such that the Carleman weight function ip 
i/>(t, x) = \x- x°\ 2 - /3t 2 + Co satisfies 

^(0,x) > Co, 35 6(0,1), 

ip{t,x)<C , te [-T,-T + ri]U[T-ri,T\, x£ (0,1). 



In particular, 



<p(P,x) > e xc °, 35 6(0,1), 

<p(t,x) < e AC °, te [-T,-T + 7/]U[T-7?,T], xG (0,1) 



In the sequel, we fix /3 as above (/? 6 (0, 1) and T^//J < sup xg ( 01 ) |x — A, so, s > such 
that Corollary [213] holds and the Carleman estimate (|2.8p holds for all h 6 (0, ho) and s 6 (so, s/h). 

Step 3. Extension and truncation. We now extend the problem (|3. 14[) on (— T, T), setting 
Zh(t) = Zh{— t) for all t G (— T, 0). We also extend dtRh in an even way and keep the same 
notations for the new problem. 

Let us define the cut-off function \ £ C°°(R; [0, 1]) such that: 



X (±T) = dt X (±T) = 
X(t) = l for all t 6 [-T + r?,T- rj\. 

We set Wh = X z h that satisfies the following equation: 

Lh[qh]w h = dttX z h + 2dtxdtz h + X fhdtRh, t G (-T, T), 

w ,h(t) = WN+i,h(t) = 0, te(-T,T), 
w h (0) = 0, dtw h (0) = f h Rh{0), 
w h (±T) = 0, d t w h (±T) = 0. 



(3.17) 



(3.18) 



Step 4. Using the Carleman estimate. From now on, C > will correspond to a generic 
constant depending on so, \,T,x° , (5, x and 77 but independent of h G (0,/io) and s G (so,e/h). 
We use the same notations as in Section [2] and set Vh = exp(sip)wh- We then have (recall (|2.39p ) 

Ph,iv h = d u v h - (1 + ylo)A h Hh + s 2 A 2 [<^ 2 (ftV 1 ) 2 - A 2 ] v h 
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and v h {±T) = d t v h (±T) = 0, v h (0) = and d t v h (0) = f h R h (0)e s ^ ^. 

Using the properties of Vh, Lemma [2 . 1 2 1 and (|2.38[) . we can make the following calculation: 

Ph,iVh d t Vh 

T J '(0,1) 

/ / {dttv h - (1 + A Q )A h v h + s 2 X 2 [ip 2 (d t ip) 2 - As] v h ) d t v h 

J -T J (0.1) 



~ 2 



|0tUh(O)f + 



(0,1) 



T J [0,1) 



9+%a+((i + Ao)a^ h ) 



s 2 A 2 



J-T J (0,1) 



-T J (0,1) 



(0,1) 



T J (0,1) 



/l/o i) (^ 9t(l ^" h|2)m ^ (1 + Ao) + ^^o^^m+a^^ 



(0,1) 

1 

2 



IM e 



s 2 C 



T J (0,1) 



T J [0,1) 



|a+^| 2 m+(a t A ) + 



T J [0,1) 



9^ Ao9^" v h m£d t v h 



J(0.1) J-' 



'(0,1) 

-Ox(sh) 



T J (0,1) 



d>TM + 



T J [0,1) 



|9tM 



T J (0,1) 



Therefore 



IM e 



2 2s<^(0,-) 



(0,1) 



T J (0,1) 

+Oa(s/i) 



7 

-T J {0,1) 



V h 



-T J (0,1) 



T J (0,1) 



Using 



Ph,iVh d t Vh 



T J (0,1) 



- ( 



|-F%,i«a| + s 



T J (0,1) 



|3 t Uh| 



T J (0,1) 



and the fact that Ox(sh) is bounded by some constant independent of s since sh < e, we get 



IM e 



2 28^(0,-) 



(0,1) 



|-F\,lM +S 



T J (0,1) 
T 



\d t v h \ 



T J (0,1) 



+ CS 5 ' 2 [ [ \v h \ 2 + CV~s([ [ \d+v h \ 2 + 
J -T J (0,1) \J-tJ[o,i) 



\d t Vh\ 



T J (0,1) 
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From the Carleman estimate ()2.43[) of Lemma fe.lSl this implies that for all s satisfying sq < s < j-, 



T r rT 



| /h |2 e 2 SV (0,) < Ml / \P h V h \ 2 + Msj \(d^V h ) N + 1 \ 

(0,1) J -T J (0,1) 



+Ms f f \hd+d t v h \ 2 

J-T J [0,1) 

< M f T [ e 2s *\L h w h \ 2 + Ms [ T e 2 ^ (M) |(9,> h ) 

J -T J (0,1) J-T 



JV+1 1 

(0,1) J-T 

+Ms f T f e 2sv \hdtd tWh \ 2 , (3.19) 

J -T J[0,1) 

where the last estimate follows from (|2.54p - (|2.55[) . 

From equation (|3. 18f> . the properties (|3.17|) of the cut-off function x an d (|3.16|) of the weight 
function ip, which is a decaying function of \t\, one gets 

[ T I e 2sv \L h w h \ 2 <C f T [ e 2 ^ {\ X ,fhd t R h \ 2 + \d tX d t z h \ 2 + \d ttX z h \ 2 ) 

J —T J (0,1) J -T J (0,1) 

< C [ T f e 2s *\f h \ 2 \d t R h \ 2 + c( f T+ \ f T ) [ e 2sv> (\d t z h \ 2 + \z h \ 2 ) 

J -T J (0,1) \J -T JT-T)/ J (0,1) 

< CK 2 f e 23V ^\M 2 +Ce 2seXa ° ( [~ T+ \ C W(t). 

J(0,1) \J-T JT-riJ 

Using now the energy estimate (|3. 15[) . 



[ e 2 ^\L h w h \ 2 < CK 2 [ e 2 ^\f h \ 2 + CK 2 e 2seXC ° [ \f h f 

•J(0.1) J (0.1) J (0.1) 



J I 1 1 

TJ(0,1) J (0,l) J (0,1) 

< CK 2 [ e 2sv{0 ^\f h \ 2 . (3.20) 

J '(0,1) 

Similarly, since dtWh = x^tZh + dtx z h, using the energy estimate (|3.15|) . 
/ e 2sv \hd+d t w h \ 2 

T J [0,1) 



2s v^ 2 \ua+a. i, i 2 _i_ oh 2 / / „ 2, via.vi 2 ifl+,. i 2 



< 2/ / e^ X Z \hdtd t z h \ 2 + 2h 2 / e 2stp \d tX \ z \dtz h \ 

-TJ[0,1) J-TJ[0,1) 



< 



2 ( T f e 2sip X 2 \hd+d t z h \ 2 + 2h 2 CK 2 f e 2 *^ \f h \ 2 , (3.21) 

J -T J [0,1) J (0,1) 



where we have used that, as proved above 

rT r / , — T+ri pT 



f f e 2 ^\d tX \ 2 \dtz h \ 2 < Ce 2 ^ C ° ( f~ + \ f ) EUt) < CK 2 f e 2s ^\f h 

J -T J [0,1) \J-T Jt- v J J (0,1) 

Therefore, plugging (pr27Jj) - (ET2Tj) in $515$ we obtain 

V~sr 2 f e 2 ^\f h \ 2 < CK 2 f e 2 ^\f h \ 2 + Cs f e 2 ^ X 2 \(d^ z h ) N+1 \ 2 

J (0,1) J(0,1) J-T 



'(0,1) J (0,1) 

+ Cs [ [ e 2sip \hd+d t z h \ 2 + Csh 2 K 2 

J-T J [0,1) -'(0,1) 



,2*V»(0.0| f. 2 
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Thus, since sh 2 < e(h)ho < 1, taking s* > so such that for all s > s„, y^sr 2 — 2CK 2 > 0, for all 
h G (0, /it) with /i» = minjfto, e/s*}, we obtain 



/ e 2s ^\f h \ 2 < Cs. f 1 e 2 «*^ X 2 \(d^z h U +1 \ 2 + Cs* f f e 2 ^\hd+d t z h \ 2 , 

J (0,1) J-T J -T J [0,1) 

and therefore 

ll/h|| L 2 (01) < C \\(d^ z h ) N+1 \\ L2( _ T T) + C \\ h dh d tz h \\ L2( _ T T;L 2^ 01)) , (3.22) 
which coincides with (|3.7[) . The proof of Theorem 13.21 is then complete. □ 



Remark 3.4. With the notations of Theorem[3~R if R h G H 1 (0, T; Lf?(0, l))nW /2,1 (0, T; L 2 (0, 1)) 
and Rh(0, ■), dtRh(0, •) G L^°(0, 1), considering the equation satisfied by Wh — dtz^: 

( dttw h - A h w h + q h w h = fhdttRh, t G (0, T), j G [1, TV], 

I w , h (t) = wjv + i, h (t) = 0, t G (0, T), 

1 to h (6) = f h R h (6), dtw h (0) =f h dtR h (Q), 

from Lemma \3.'A we get 

sup E%(t) < CE)"(Q) + C \\fh\\ 2 L 2 (0 1} ||-Rh||^2,i( t-l 2 (o D) ' 
te(o,T) h - J y hK " 



with 



Et(Q)=f \d+(f h R h (0))\ 2 + [ \f h d t R h (0)\ 2 + [ \f h R h (0)\ 2 

J\0,1) J (0,1) ■'(0,1) 

d+(f h )\ 2 \m+R h (0)\ 2 + C \\f h \\l l(0A) \\(R h (0, ■), d t R h (0, -))||i- co ,i) a 



< 2 



'[0,1) 7(0,1) 7(0,1) 

I";. U/. 'I "/>'")'" 

'[0,1) 

Therefore, if in addition to (|3. 4[) . there exists a constant S > smc/i i/wrf for all h > 0, 

II-RaIIw^. 1(0, T;ij(0,l)) + II (^Rft (0, C?t-R/i (0' ')) llz,?^ (0,1)2 < ~r , 

V r (3-23) 

\dih\ 2 <^ im 2 , 

[0,1) " 7(0,1) 

tften 

W A\ L ~(0,T;Ll[0,l)) = \\ d h 9 t Z h\\L°°(Q,T;Ll[0,l)) ^ C5fl ~ 1 \\Ml%(P,1) ' 

and in particular, 

h\\d h dtZh\\ LOO(!0tT . L a [0>1)) < C<5 \\fh\\ L a (0,1) ■ 

Therefore, if condition (|3.23|) is satisfied for S > small enough, estimate (|3.22|l simply becomes, 
for h small enough, 

ll/h|| L 2(0,l) ^ C IK^h Z ' l ) Ar + 1 llL2 ( _ TjT) • 

Condition (I3.23|l can fee seen as a filtering condition on the data. To be more precise, if we filter 
enough the data (at the scale 8/h with 5 small enough ), the Tychonoff regularization term is not 
needed anymore in (|3.7p . 

3.3 Uniform stability for the discrete inverse problem 

Proof of Theorem \3.1\ Setting uh = Vh[qh] — Vh\ph], where yu[qh] and yh[ph] are respectively the 
solutions of (|1.7[) corresponding to p^ and g^, then Uh solves 

f dttu h - A h u h + q h u h = f h R h , t G (0, T), j G [1, AT], 

I u ,h(t) = u N +i,h(t) = 0, t G (0, T), (3.24) 

{ u h (0) = 0, d t u h (0) = 0, 

with fh=Ph — qu and Rh = yh\ph]- We then directly apply Theorem 13.21 □ 
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Remark 3.5. Remark \3-4\ also applies here of course, and the filtering condition (|3.23p then 
becomes: 

bh[Ph]lll V 2,l(0,T:L2(o.l)) + l|(yh^h)|ll,-(01)2 ^ T> / \ d h (ih - Ph)\ 2 < / \lh ~ Phf , 

ft h K ' ' n J[o,i) n J(o,i) 

for some 5 > small enough. A convenient way to satisfy these two conditions is to impose that 
both ph and qu belong to a filtered space and that the data (j/h,J/h), gn and (<?h,</h) are smooth. 

4 Convergence issues 

In this section, we will detail and prove the convergence results that were presented rapidly in 
the introduction. 

4.1 Statements of the results 

In order to prove a convergence result, we shall need some assumptions first. 

Assumption 1 (A priori bounds on the potential). There exists m > such that p G L^ m (0, 1). 

Assumption 2 (Regularity assumptions). The data satisfy 

O/V^tf^O^xff 1 (0,1), 
g G W 1 ' 1 (0,T;L 2 (0,1)), (ffV) G (H 2 (0, T)) 2 , 
with the compatibility conditions 

g°(0) = 2/ (0), ^(0) = y°(l), dtg°(0) = t/(0) and 8^(0) = j/(l). 

One should notice that under these regularity assumptions, according to [27] (see also Re- 
markEU), for p G L< m (0, 1) the solution y[p] of (fTTTJ) belongs to the space C 2 ([0, T]; L 2 (0, 1)) n 
C 1 ([0,r];fl" 1 (0,l))nC ([0,T];i7 2 (0,l)). In particular, one can check that d tx y[p](-, 1) G L 2 (0,T ) 
(this can be found in [27] but can also be seen as a consequence of the multiplier identity l|4.12[l ) 
and that y[p] G H 1 (0, T; L°°(0, 1)). 

Since we are interested in a convergence result, we shall explain how to compare discrete 
functions with continuous ones. In order to do so, we introduce two extension operators. 

The first one extends discrete functions by continuous piecewise affine functions. To be more 
precise, if fh is a discrete function (fj,h)jeio,--- .iv+lji the extension eh(fh) defined on [0, 1] by 

Ch(fh)(x) = Un + ( L ! h 7 '•) (x - jh) on [jh, (j + l)h], j G 10, ■ ■ ■ , Nj. 

This extension presents the advantage of being naturally in // 1 (0, 1). 

The second one is the piecewise constant extension e°(/h), defined for discrete functions 

{fj,h)jE{l,- ,nj by 

eUA) = fi,h on [(j - l/2)h, (j + l/2)h[, j G [1, • ■ ■ , Nj, 
el(f h ) = on [0, h/2[U[(N + l/2)h, 1]. 

Of course, this one is more natural when dealing with functions lying in L 2 (0, 1). In particular, 
we have 

||eh(/h)|| L 2 (0 , 1) = ll/fc|lz,a(o,l) • 

Also note that easy (but tedious) computations show that eh(fh) converge to / in L 2 (0, 1) if 
and only if e°(/h) converge to / in L 2 (0, 1). 

Of course, we shall need some convergence estimates: 
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Assumption 3 (Convergence assumptions). The sequence of discrete data (yh,Vh) satisfies 

(e£(A fc y£),e£(A fc tfi)) — > (Ay , Ay 1 ) L 2 (0, 1) x H-\0, 1). (4.1) 



27ie sequences of source terms gh, (ghi9h) satisfy 

e° h (gh)^9 in W M (0,T;X 2 (0,1)), (<?£, y, 1 ) — »• (<?°, ff 1 ) in (H 2 (0, T)) 2 . (4.2) 

Finally, we shall also need a uniform positivity assumption: 
Assumption 4 (Positivity). There exists r > swc/i t/iat 

inf {\y°(x)\,x 6 (0, 1)} > r > 0, and V7i > 0, ^inf^ |y°,J > r. (4.3) 
Now, we introduce, for h > 0, the following observation operator: 

h : L£< m (0,l) -> L 2 (0,T) xL 2 ((0,T) x (0,1)) 

Ph i-> (^t(9 h i/h[ph])jv+i, hd x e h (d u yh[ph])) , 

where yh\ph] is the solution of (|1.7|) with potential ph- We also introduce its continuous analogous 

0o: Lf ro (0,l) -> L 2 (0,T)xL 2 ((0,T)x(0,l)) 

P ^ (9Ay[p](-,l)>0), l4 '° j 

where y[p] is the solution of (|l.l|l with potential p. 

Note that, using these notations, Theorem 13. II can then be seen as a uniform stability of the 
maps G^ 1 . Indeed, 

\\hd£dttyh[qh] - hd hdttyh[ph}\\ L2{QT . L 2 h{[01)) 

= \\hd x e h (d t tyh[qh]) - hd x e h (dttyh[ph})\\ L 2 ((0tT)> , (0A)) , 

and then (|3.3[) reads as: 

||e°(?h) - e h(Ph)|| L 2 (0il) < C \\&h(ph) - ©h(<7ft)llL2(o,T)xL2((o,T)x(o,i)) • (4-6) 
Our main result is then the following convergence theorem: 
Theorem 4.1. Under Assumptions]!}^ let q^ £ L^ <m (0, 1) be such that 



@h(qh) — + 6o(p) strongly in L 2 (0,T) x L 2 ((0,T) x (0, 1)). (4.7) 

h— J-0 

Then one has the convergence 

e° h (q h )^P *nL 2 (0,l). (4.8) 

h— ¥0 

Before going into the proof of Theorem 14.11 we shall emphasize that there exist discrete 
sequences of potentials such that (14.711 holds. Actually, this is a consequence of the following 
consistency result: 

Theorem 4.2. Under Assumptions Qr[3 f or a ^ potential p 6 L^? m (0, 1) there exists discrete 
potentials pu £ L^ <m (0, 1) such that 

e° h (p h )^P mL 2 (0,l) and Q h (p h ) — ► 9 (p) mL 2 (0,T) x L 2 ((0,T) x (0,1)). (4.9) 

h— >0 ft— >0 

Moreover 

SUP ||yh.[pfe]||H- 1 (0,T;£°°(0,l)) < °°' ( 4 - 10 ) 

fee (o.i) 

where yh[Ph] is the solution of ((TTTJ) . 

In the following section, we give the proofs of these Theorems. Actually, as we will see, 
Theorem l4.2l is the second milestone of the proof of Theorem 14. II the first one being Theorem l3.ll 
In other words, the proof of Theorem 14. II that will be given at the end of this section, relies on a 
Lax-type argument for the convergence of the numerical schemes based on the consistency of the 
method, given by (|4.9[) . and the uniform stability (|4.6[) of the discrete inverse problems. 
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4.2 Proofs 



Proof of Theorem \4-S\ In the proof, we shall distinguish the regularity and convergence issues 
coming from the boundary source terms and the initial data from the classical ones coming from 
the potential and distributed source term. 

Step 1: Convergence without potential and source term. 

Let z be the solution of 

( d a z - d xx z = 0, (t, x) £ (0, T) x (0, 1), 

{ z(t,0) = g°(t), y(t,l)=g 1 (t), t € (0,T), (4.11) 

[ a(0,0 = w°, dtz(0,-)=y\ 

Since (i/ ,?/ 1 ) and (g ,^ 1 ) satisfy Assumption [2] the solution z of (|4.1ip satisfies (see [27] for 
details) 

dtz G C([0, T];H\Q, 1)) n C\[Q, T], L 2 (0, 1)). 
We are therefore allowed to write the following multiplier identity 

i [ T \d xt z{t,l)\ 2 dt = i / T A 1 (|«9 tt z| 2 + IcM 2 ) dxdi-J r^tg^dt 
^ Jo * Jo Jo ^ J 

+ [ d u z(T,x)xd xt z(T,x)dx- [ d u z(0,x)xd xt z{0,x)dx, (4.12) 
Jo Jo 

which is obtained by differentiating in time equation (|4.11|l . and then multiplying it by xdt x z, 
integrating over (0, T) x (0, 1) and doing integration by parts. 
Now, let Zh be the solution of 

[ duz^n - (A h z h ) 3 = 0, t6(0,T), j€{l,N\, 

I z Q . h (t)=g° h (t), z N+1M (t) = gi(t), t€[0,T\, (4.13) 

{ z jth (p) = yl h , d tZj , h (0) = Vlh, j G [1, Nl 

In this step, we want to prove that Zh — > z in the appropriate functional spaces. In order to do 
that, we use the following result: 

Theorem 4.3 ( 17 ). Let fl)h>o be a sequence of boundary data strongly convergent to some 
functions (f°,f l ) in L 2 (0,T) 2 . Let {<p° h , ip\) be a sequence of discrete functions such that 

(eU<Ph)A(<Ph))r+ n (<(' >'f >1 ) st ™ngly m L 2 (0, 1) x H' 1 (0, 1). (4.14) 
Then the solutions Lph of 

f du(Pi,h - (A^) 3 . = o, te (o,T), j e li, N}, 

{ **>,*(*) = .#(*), <PN+i,h(t) = fk(t), * 6 [0,7], 

[ w , h (o) = ftw.hCo) = v) th , j g [i, iv] 

converge toward the solution ip of 

( dtt<p - d xx ip = 0, (t, a) G (0, T) x (0, 1), 

I p(t,0) = /°(t), ^(t,l) = / 1 («), te(o,T), 
I ^(0,-) = ¥>°, &p(M = <p\ 

in tfte following sense: for all p < oo, 

e°(y>h) — > V ^ron^y l nL p ((0,r);L 2 (0,l))nH/ 1 ' p ((0,T); J ff- 1 (0,l)). (4.15) 

h— >0 

Besides, for all to G [0, T], 

(e° h (<p h ){t ),d t e h {if h )(to)) — ► (<p(to),dt<p(to)) strongly in L 2 (0, 1) x i2 _1 (0, 1). (4.16) 
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For the proof of Theorem 14. 31 we refer to [17]. Note that Theorem 231 is not standard, since it 
deals with solutions of the continuous wave equation defined in the transposition sense. Therefore, 
the proof of Theorem 14.31 is based on a duality argument and convergence results for the adjoint 
equation, namely the waves, and in particular on their normal derivatives on the boundary (which 
corresponds to the adjoint operator of the Dirichlet boundary conditions). 

Of course, regarding the regularity hypothesis in Assumption [2] and the convergence one in 
Assumption [3] we can apply this result to z^, dtZh and dttZh- Of course, the latter yields the 
strongest result, thus improving the ones on dtZu and z^: for all p < oo, 



f dtte° h (z h ) — » dttz strongly in L p ((0, T); L 2 (0, 1)) n ^((0, T); #^(0, 1)), 

d t e h (z h ) ^> 8 t z strongly in L p ((0,T); J ff 1 (0,l))nW 1 ' p ((0,T);L 2 (0,l)), ( 4 .17) 
e h {z h ) z strongly in W 1 '" ({0, T); H 1 (0,1)) f]W 2 ' p {(0,T); L 2 {0,1)), 



and, for all t € [0,T], 

(e o h (d h d t z h )(t ),e° h (dttz h )(t o )) — ► (d xt z(t ),dttz(t )) strongly in (L 2 (0, l)) 2 . (4.18) 

h— >0 

Now, we focus on the convergence of the normal derivatives. This is slightly more subtle. 
First, arguing as in [T7J by duality against smooth functions, one easily checks that 

(d^z h ) N +i d x z(;l) weakly in L 2 (0,T). (4.19) 

h— >0 

Then, we derive a multiplier identity similar to (|4.12[) for the discrete equation (|4.13|) . In order 
to do that, we multiply equation (|4.13l) differentiated once in time by xdhdtZh' 

I dtttz h x d h d t z h dt- / A h d t z h x d h d t z h dt = 0. 
o J(o,i) Jo J(0,1) 



But, on the one hand, using (|2.19[) . we have 



/ dtttz h x d h d t z h dt = d tt Zh x d h d t z h 

J(0,1) J(0,1) 



T T 

— / dttz h x dhdttzn dt 
o Jo J (0,1) 



dttz h x d h d t z h 

(0,1) 



o 

and on the other hand, using now (|2.22[l . we get 



T 

+ \ [ [ \duz h \ 2 dt-£ f f \d+d tt z h \ 2 dt 

z Jo J (0,1) 4 Jo J [0,1) 



f A h d t z h x d h d t z h = -- / \d£d t z h \ 2 + -\(d h d t z h ) N+1 \ 2 . 

J(0,1) Z J [0,1) 



'(0,1) A J[o,i) 
Combining these last three identities, we obtain 
rT j,2 fT 

2 



\ r \(d^d t z h ) N+1 \ 2 dt+^ r f \d+d tt z h f 

z Jo 4 Jo J [0,1) 



dttZh x dhdtz h 

(0,1) 



T 

T 







+ \ ! \dttz h \ 2 dt+\ I I \d+d t z h \ 2 dt. (4.20) 
1 Jo J (0,1) 1 Jo J[0,V) 



According to the strong convergences in (|4.17[) and (|4.18[) , we can pass to the limit in the right 
hand side of (|4.20[) , which converges to the right hand side of (|4.12|) . leading to: 

lim (i / \(d^d tZh ) N+1 \ 2 dt + ^j- f I \d+d tt z h \ 2 dt\ =i / \d xt z(-,i)\ 2 dt 



hTo \ 2 L .v-« /--r-i - 4 



o J [0,1) 
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This last fact, together with the weak convergence (|4.19|l . implies that 
((d^dtz h ) N+1 ,hd x e h (dttz h )) — ■+ (d xt z(; 1), 0) strongly in L 2 (0,T) x L 2 ((0,T) x (0,1)). (4.21) 



Step 2 Convergence with source term and potential 

So far, we did not assume anything on the potential p and on the convergence of the discrete 
potentials ph to p, since they did not appear in the study of z and Zh solutions of (|4.11[) and (|4.13|) . 

Since p £ L^ m (0, 1), it is very easy to construct a sequence ph £ L^ <m (0, 1) such that e1(ph) 
strongly converge to p in L 2 (0, 1). Taking such sequence ph, we set yh[Ph] = Zh + Vh[ph] where 
^h[ph] is the solution of 



Due to the convergence hypothesis in Assumption [3] we have the convergence of e ( l(gh) towards 
g in W 1 ' 1 (Q,T;L 2 (0, 1)) as h 0. 

From (|4.17[1 . e°(zh) and e1(d t Zh) strongly converge in L 2 (0, T; L 2 (0, 1)) towards z and 9tZ, 
respectively. Besides, eh{zh) and eh{dtZh) respectively converge to z and dtz in L 2 (0, T; i/ 1 (0, 1)). 
Since z ,fc(i) = flh(t) and d t z ,h{t) = d t ff°(t) are bounded (# 2 (0,T) C C^QO.T])) and strongly 
converge in C°([0,T]), respectively, toward z(i,0) = g°(t) and <9tz(t,0) = dtg°(t), eh{zh) and 
&h{dtZh) respectively converge to z and dtz in L 2 (0, T; L°°(0, 1)). Of course, this implies the 
L 2 (0, 7 1 ; L°°(0, 1)) boundedness of the sequences e° and e° (dtZh)- Since they converge strongly 
in L 2 (0, T; L 2 (0, 1)), we also have the strong L 2 (0,T; L 4 (0, l))-convergences of (e < l(zh),e < l(dtZh)) 
toward (z, d t z). 

Similarly, e°(p k ) strongly converges to p in L 2 (0, 1) and is bounded in L oo (0, 1). Therefore, 
e h(Ph) strongly converges to p in L 4 (0, 1). 

Since eKanbh) = e° (afc)e°(6^), we thus obtain that e1(phZh) and eKphdtZh), respectively, 
strongly converge to pz and pdtz in L 2 (0, T; L 2 (0, 1)). 

Therefore, 

el(g h -p h z h )—>g-pz in W /1 ' 1 (0, T; L 2 (0, 1)). 

h— >0 

Thus, classical results yield the convergence of Uh[ph] toward v[p], solution of 
( dttv- Av+pv = g~pz, (t,x) e (0,T) x (0,1), 



Therefore, we obtain 

J fte k (» fc |p fc ])^ft»|p] strongly in L 2 ((0,T); ^(0,1)) n /^((O, T); L 2 (0, 1)), 

I e fc (»fclpfc]) — Ti] strongly in L 2 ((0,T); ^(0,1)) n^ 2 ((0,r);L 2 (0,l)) (4 ' 22) 

and, for all t G [0,T], 

(e°(9 k at« k [p fc ])(to),e°(9tt?; fe [pfe])(to)) — > (d xt v[p](to),dttv[p](t )) strongly in (L 2 (0, l)) 2 . 

h >l! 

Of course, as for Zh, using the discrete multiplier identity satisfied by dtv (see (|4.20p ^ and the 
above convergences, we easily get 

((dhd t v h [p h ]) N+1 ,hd x e h (d u v h [p h ])) — ► (d xt v[p](- , 1) , 0) 

h->0 

strongly in L 2 (0, T) x L 2 ((0, T) x (0, 1)). (4.23) 

Now, using (|4.21l) and (|4.23l) . the solution yh[ph] of (|1.7|) converges to the solution y[p] of 
in the following sense: 

((d^d t y h [p h ]) N+1 ,hd x e h (duy h [p h ])) — > (^ t y[p](-,l),0) strongly in L 2 (0, T) x L 2 ((0, T) x (0, 1)), 





34 



which is precisely (14. 9p . 

Besides, using (14. 1 Tf) and (14.221) . we have 

e h (y h \p h ]) — > y[ P ] in ff^O.Tjff^O.l)), 

which of course implies the bound (|4.10|) . This concludes the proof of Theorem 14.21 □ 

We are now in position to prove Theorem 14. II 

Proof of Theorem\Jl\ Let p € L< m (0, 1) and let q h € L~< m (0, 1) be such that ((4~7)) holds. 
Denote by ph the potentials given by Theorem 14.21 Then we have 

Q h {Ph) - <dh(q h ) — > (0,0) strongly in L 2 (0,T) x L 2 ((0,T) x (0,1)). 

But according to (|4.10p and the positivity Assumption 2] we can apply Theorem 13.11 for some 
C > independent of h > 0, estimate (|4.6p holds. Therefore, e°(p rt ) — e^(qh) strongly converges 
to zero in L 2 (0, 1). Using (|4.9[) . we deduce that e°(gh) strongly converges to p in I/ 2 (0, 1). □ 

5 Further comments 

• Other convergence results. Note that our convergence results require the convergence 
of hd^dttVh[Ph] to zero in L 2 (0, T; L 2 (0, 1)). This term is here to handle spurious high-frequency 
waves generated by the space semi-discretization - see e.g. [32] - which are by now well-known 
to be responsible for the lack of uniform observability of waves [36]. Of course, other ways of 
removing these high-frequency waves can be implemented, an easy one being to impose some 
smoothness and filtering conditions on the data - see Remark 13.51 Note however that these 
conditions seem to be more difficult to implement in practice. 

• Time discretization. Here we focused on the space semi-discretization of the wave equa- 
tion for simplicity. Indeed, the fully discrete wave equation in which the time-derivative has been 
approximated by the centered difference approximation could be handled the same way, since 
time and space are completely decoupled then. This will of course introduce a Tychonoff regu- 
larization term within the Carleman estimates of the same order but depending not only on the 
space discretization parameter, but also on the time semi-discretization parameter. This again is 
completely compatible with the known results on the observability of discrete waves - see [151117] . 

• Other space discretizations. Here, we have chosen a very simple space discretization 
process corresponding to the finite-differences approach. Other space discretizations should be 
studied, but regarding the literature in what concerns discrete observability estimates for the 
waves (see e.g. [361 116] ). we expect the Tychonoff regularization term to be needed within the 
discrete Carleman estimates in the case of finite-elements methods. However, for mixed finite 
elements methods (see [101 1111 113]). we expect better behavior than here and this Tychonoff 
regularization term may be not needed anymore. This should be studied carefully. 

• Higher dimensions and more sophisticated wave models. Of course, an interesting 
question would be to develop these discrete Carleman estimates in higher dimensions (as it was 
done in [7])- It is usually admitted that Carleman estimates "do not see" the dimension of the 
space. This is indeed true in the continuous case, but in the discrete case, the integrations by 
parts are much more intricate. This is currently under investigation. Regarding more generic 
hyperbolic models, one could also mention [24], [4] or [2] giving stability of inverse problem from 
global Carleman estimates respectively for the Lame system, a discontinuous wave equation or in 
a network of 1-d strings. 

• Semilinear wave equations. One of the standard applications of Carleman estimates 
is to prove controllability of semilinear wave equations - see |12l 118] . We expect that these 
discrete Carleman estimates could be of some use to prove the convergence of discrete controls 
for semilinear wave equations and to improve the results already obtained for globally Lipschitz 
nonlinearities using bi-grids methods in [37] . 
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• How to compute a discrete sequence ph such that @h(ph) converges to ©o(p) ? 
This is certainly one of the most challenging issues concerning this kind of inverse problems, since 
the map Oh is highly nonlinear. Of course, a natural idea is to introduce 

Jh(Ph) = \\@h(ph) ~ &o(p)\\l3(p,T)xl'(,o,t ; l'(p,i)) 

and to minimize it on the set L£° <m (0, 1). But this can be very hard since Jh may have several 
local minima. Another approach will be presented in the work [3] based on Carleman estimates 
and stability results inspired from [221 IT] . 

Acknowledgements. The authors acknowledge Jerome Le Rousseau, Franck Boyer, Jean-Pierre 
Puel, Masahiro Yamamoto and Frederic De Gournay for interesting discussions related to that 
work. 
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